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EXECUTIVE  SUMMARY 


This  program  was  initiated  to  investigate  the  potential  of 
nuclear  magnetic  resonance  (NMR)  for  remote  detection.  Remote 
regions  accessible  from  only  one  side  and  the  remote  detection 
"  of  explosives  were  topics  of  particular  interest. 

The  program  involved  the  following  seven  specific  tasks  and 
^  subtasks:  (1)  to  identify  the  compounds  of  interest,  both  benign 

and  explosive;  (2)  to  locate  and  catalog  available  data  on  the 
compounds  of  interest;  (3)  to  review  current  remote  NMR  technol¬ 
ogy;  (4)  to  identify  alternative  remote  NMR  systems;  (5)  to 
define  a  f igure-of-merit  suitable  for  comparing  the  alternative 
NMR  systems;  (6)  to  use  the  f igure-of-merit  to  identify  the  most 
promising  of  the  alternative  systems;  and  (7)  to  predict  the 
performance  of  the  chosen  system.  Each  of  these  tasks  was  com¬ 
pleted  successfully. 

The  present  state  of  remote  NMR  technology  and  the  principal 
conclusions  of  this  study  are  outlined  below. 


*  Remote  detection  by  single-sided  NMR  is  feasible. 
Present  systems  produce  low  static  fields  (less  than 
500  Gauss)  at  modest  penetration  depths  (10  cm). 

*  The  most  promising  advanced  remote  NMR  system  con¬ 
figuration  utilizes  an  Inside-Out  Helmholtz  pair  of 
current  loops  to  produce  the  static  magnetic  field 
and  a  semitoroid  coil  to  produce  the  rf  field. 

*  The  advanced  NMR  remote  detection  system  will 
employ  superconducting  coils  to  produce  magnetic 
fields  of  1,000  Gauss  or  more  at  penetration  depths 
in  excess  of  20  cm.  New  cryogenic  materials  and 
fabrication  techniques  should  eliminate  many  of  the 
difficulties  previously  associated  with  the  use  of 
superconducting  technology. 

*  Available  data  indicate  that  NMR  remote  detection 

of  the  following  explosives  is  possible:  Compound  B, 
Compound  C-4,  PETN,  RDX,  TAT3,  and  TNT.  Data  on  the 
NMR  properties  of  explosives  are  sparse. 

*  It  may  be  possible  to  define  a  unique  signature 
for  many  explosive  compounds  by  taking  advantage  of 
hydrogen-nitrogen  nuclear  interactions. 


Additional  laboratory  data  and  analyses  are  required  to  determine 
the  range  of  applicability  of  this  technology  and  to  address  the 
practical  proble.ms  associated  with  its  implementation. 
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NTRODUCTION 


I . 


The  reliable  detection  of  explosives  in  remote  or  inaccessible 
regions  poses  difficult  technical  problems.  The  purpose  of  the 
present  program  is  to  explore  and  evaluate  the  potential  of  one¬ 
sided,  or  single-sided,  nuclear  magnetic  resonance  ( NMR )  for  the 
remote  detection  of  explosiv  s. 

NMR  utilizes  nonionizing  radiation  to  manipulate  the  weak  magnet¬ 
ism  of  atomic  nuclei.  The  non-ionizing  radiation,  a  combination 
of  AC  and  DC  magnetic  fields,  is  applied  to  the  region  to  be 
inspected.  Shortly  thereafter,  a  radio  frequency  signal  is 
received  from  the  nuclei.  The  characteristics  of  this  signal  are 
related  to  the  chemical  properties  of  the  material  containing  the 
nuclei . 

This  study  was  limited  exclusively  to  NMR  of  hydrogen  nuclei. 
Hydrogen  nuclei  give  a  strong  NMR  response  and  are  abundant  in 
most  of  the  materials  of  interest.  Hydrogen  is,  without  doubt, 
t.he  only  nucleus  likely  to  provide  an  NMR  signal  strong  enough 
for  the  detection  applications  considered. 

In  addition,  under  certain  conditions,  the  hydrogen  NMR  signals 
of  many  explosive  compounds  possess  an  al.most-unique  signature. 
This  signature,  based  upon  differences  in  NMR  parameters  called 
relaxation  times,  provides  the  foundation  for  a  simple  and  direct 
detection  method. 

NMR  detection  systems  are  capable  of  performing  well  in  a  field 
environment.  As  an  example,  an  NMR  system  to  detect  narcotics  in 
the  mails  was  recently  field  tested  at  the  Worldway  Postal  Center 
in  Los  Angeles  (1).  The  NMR  system  was  designed  to  respond  to 
the  relaxation  signature  of  cocaine.  A  total  of  3,350  letters 
were  examined.  Seven  alarm  responses  were  recorded;  two  were 
false  alarms,  and  five  were  due  to  letters  containing  cocaine. 

An  important  reason  for  the  success  of  this  field  test  was  that 
the  NMR  system  was  not  forced  to  encounter  large  quantities  of 
benign  and  potentially-interfering  materials  (other  t.han  paper) 
during  the  inspection  process.  Since  t.he  volume  of  the  inspec¬ 
tion  region  in  remote  detection  NMR  is  limited,  it  should  be 
noted  that  the  NMR  systems  discussed  in  this  report  may  display 
similarly  favorable  inspection  properties,  even  in  complex  de¬ 
tection  situations. 

T.his  study  extends  previous  work  in  remote  detection  NMR, 
principally  Single-Sided  NMR  and  Inside-Out  NMR.  Previous 
i.mplementations  of  Single-Sided  NMR  (2,3,4)  have  employed  iron- 
yoke  electromagnets  to  produce  the  required  region  of  re.mote 
field.  These  systems  were  developed  to  measure  moisture  levels 
in  subsurface  soils  and  in  structural  concrete. 
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Inside-Out  NMR  (5,6,7)  utilizes  opposed  magnets  to  produce  a 
toroidal  field  geometry.  This  type  of  remote  NMR  system  was 
developed  for  geophysical  exploration,  i.e.,  as  a  borehole 
logging  tool  to  evaluate  subsurface  rock  formations  and  hydro¬ 
carbon  deposits. 

The  principal  objectives  of  this  study  were  to  define  the  problem 
of  NMR  remote  detection  of  explosives;  gather  available  data  on 
the  materials  of  interest,  both  benign  and  explosive;  define  a 
figure-of -merit  suitable  for  comparing  alternative  instrument 
configurations;  and  utilize  the  f igure-of-merit  to  identify  one 
or  more  remote  detection  NMR  systems  with  the  potential  to  con¬ 
tribute  to  the  solution  of  the  detection  problem.  All  of  the 
principal  objectives  of  this  study  were  achieved. 

This  report  is  organized  along  the  lines  of  the  study  itself, 
beginning  in  Chapter  II  with  a  brief  presentation  of  NMR 
fundamentals.  A  discussion  of  the  NMR  properties  and  molecular 
structures  of  the  materials  of  interest,  which  include  both 
explosive  and  benign  substances,  follows  in  Chapter  III. 

A  discussion  of  system  components  for  remote  NMR  is  presented  in 
Chapter  IV.  A  review  of  the  possible  combinations  of  these 
components  leads  quite  naturally  to  several  potential  systems. 
These  systems  are  discussed  and  evaluated  in  Chapter  V. 

The  most  promising  configuration  for  a  single-sided  NMR  system  is 
presented  in  Chapter  VI.  The  system  signal-to-noise  ratio  is 
calculated  and  practical  considerations  are  discussed.  Conclu¬ 
sions  and  recommendations,  which  are  quite  positive  concerning 
the  potential  of  single-sided  NMR  for  explosives  detection,  are 
presented  in  Chapter  VII. 
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II.  NMR  FUNDAMENTALS 


A.  Introduction 

The  phenomenon  of  nuclear  magnetic  resonance  (NMR)  arises  from 
the  fact  that  the  nuclei  of  many  atoms  possess  both  magnetic 
moments  and  angular  momenta  (8,9).  In  particular,  hydrogen 
nuclei,  or  protons,  which  are  found  in  most  materials,  exhibit 
magnetic  properties  especially  suitable  for  the  observation  of 
nuclear  magnetic  resonance. 

When  subjected  to  a  static  magnetic  field,  the  hydrogen  nuclei 
align  at  an  angle  to  the  field,  pointing  either  with  the  field  or 
against  it,  and  then  precess  about  the  field  at  a  frequency 
called  the  Larmor  frequency.  The  static  magnetic  field  is  often 
called  the  Ho  field.  The  Larmor  frequency  is  proportional  to  the 
strength  of  the  applied  field,  and  for  an  applied  field  of  1,000 
Gauss  (about  2,000  times  the  strength  of  the  eart.h's  magnetic 
field)  the  Larmor  frequency  is  4.26  Megahertz  (MHz).  For 
readi ly-attainable  magnetic  fields,  and  for  most  nuclei,  Larmor 
frequencies  fall  within  the  radio-frequency  (or  rf )  range  from 
about  0.3  -  300  MHz. 

Nuclear  magnetic  resonance  arises  when  an  oscillating,  or  rf, 
magnetic  field  is  applied  to  the  precessing  nuclei.  This  rf 
field  is  called  the  Hi  field,  and  it  is  applied  in  a  direction 
perpendicular  to  the  Ho  field. 

If  the  Hi  field  oscillates  at  the  Larmor  frequency,  then  some  of 
the  nuclei  that  were  originally  aligned  with  the  static  field 
will  realign  against  it.  The  net  effect  is  the  resonant  absorp¬ 
tion  of  a  very  small  cunount  of  energy  from  the  rf  field,  and  this 
resonant  absorption  is  called  nuclear  magnetic  resonance  (NMR). 


B.  Pulse  NMR 

A  common  way  to  observe  a  nuclear  magnetic  resonance  signal 
involves  placing  a  sample,  which  contains  the  nuclei  of  interest, 
into  a  static  magnetic  field  and  then  applying  a  short,  intense 
pulse  of  rf  field  at  the  Larmor  frequency.  The  technique  is 
called  pulse  NMR  (10). 

Following  the  rf  pulse,  the  nuclei  precess  in  phase.  If  a  coil 
is  located  nearby,  this  in-phase  precession  of  nuclear  magnets 
induces  a  transient  voltage  across  it.  This  transient  voltage  is 
the  NMR  signal,  which  is  cdlled  the  free  induction  decay  (FID). 

A  sketch  of  an  NiMR  FID  is  shown  in  Figure  II-l.  The  adjacent 
coil,  often  called  the  sa.mple  coil,  is  generally  wound  around  the 
sa.mple  to  ensure  efficient  coupli.ng.  The  sample  coil  is  also 
used,  in  many  configurations,  to  apply  the  rf  pulse. 


Figure  IT- 


The  magnitude  of  the  FID  signal  following  the  rf  pulse  depends 
upon  the  duration  of  the  pulse.  A  pulse  of  a  certain  length, 
called  a  90-decree  pulse,  produces  a  maximum  in  the  transient 
in-phase  precession  of  the  nuclear  moments  and,  hence,  a  maximum 
FID  signal  following  the  pulse.  In  fact,  the  effect  of  a  90- 
degree  pulse  is  to  tip  the  nuclear  magnetization  away  from  the 
direction  of  the  Ko  field  by  90  degrees. 

A  pulse  of  twice  the  length,  called  a  180-degree  pulse,  produces 
a  minimum  in  the  transient  nuclear  precession  and,  under  ideal 
conditions,  results  in  a  zero  FID  response.  In  a  fashion  similar 
to  the  90-degree  pulse,  the  effect  of  a  180-degree  pulse  is  to 
tip  the  nuclear  magnetization  away  from  the  direction  of  the  Ho 
field  by  180  degrees. 

Short,  intense  rf  fields  impressed  upon  the  nuclear  spin  system 
produce  effects  that  are  difficult  to  characterize  in  a  brief 
discussion.  For  an  in-depth  treatment  of  the  subject,  compre¬ 
hensive  references  are  readily  available  (10,11). 


C.  Laboratory  Pulse  NMR  Systems 

A  block  diagram  of  a  typical,  laboratory-type  pulse  NMR  system  is 
shown  in  Figure  II-2.  The  magnet  applies  a  static  magnetic  field 
to  the  sample,  while  the  radio  frequency  transmitter  generates 
the  rf  pulses.  The  rf  pulses  are  applied  to  the  sample  through 
the  coupling  network  and  sample  coil,  which  are  also  used  to 
receive  the  NMR  free  induction  decay  (FID)  signal. 

The  NMR  receiver,  tuned  to  the  Larmor  frequency,  amplifies  the 
weak  FID  signal  prior  to  processing  and  display.  The  controller 
is  responsible  for  circuit  timing?  it  controls  t.he  duration  of 
the  pulses,  the  number  of  pulses  in  the  test  sequence,  the  timing 
of  data  acquisition  windows,  and  other  similar  tasks. 


D.  Remote  Pulse  NMR  Systems 

The  block  diagram  of  an  NMR  system  for  remote  detection  is  very 
similar  to  that  of  the  laboratory  system  described  above.  How¬ 
ever,  since  access  to  the  sample  is  restricted,  remote  NMR 
i.mposes  special  requirements  on 'both  the  magnet  system,  and  the 
sample  coil. 

There  are  t.nree  problems  in  implementing  single-sided  remote  NMR: 
(1)  a  magnet  geometry  must  be  defined  that  will  produce  a  homo¬ 
geneous  DC  field  in  a  re.mote  region  where  access  is  limited  to 
only  one  side;  (2)  an  rf  field  must  be  produced  within  tnis 
remote  region  of  DC  field;  and  (3)  the  direction  of  the  rf  field 
must  be  perpendicular  to  the  direction  of  the  DC  field.  Practi¬ 
cal  solutions  to  each  of  three  problems  will  result  in  a  useable 
re.mote  detection  NMR  system. 


Block  Diagram  of  a 
Pulse  NMR  System 


Fi<?ure  II-2 
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E.  Spin-Lattice  Relaxation 

In  pulse  NtMR,  if  the  rf  pulse  is  of  the  proper  duration  and 
intensity,  the  nuclei  absorb  energy  during  the  pulse  and, 
therefore,  at  the  end  of  the  pulse,  are  out  of  equililbrium  with 
their  surroundings.  The  characteristic  time  for  the  nuclei  to 
return  to  equilibrium  is  called  the  spin-lattice  relaxation  time, 
Tl.  For  protons  in  common  substances  in  typical  magnetic  fields, 
Tl  ranges  from  hundreds  of  microseconds  to  hundreds  of  seconds. 

Spin-lattice  relaxation  times  are  commonly  measured  in  the  lab¬ 
oratory  with  the  (130  -  tau  -  90)  pulse  sequence  shown  in  Figure 
II-3.  First,  a  180-degree,  or  "pi'*,  pulse  is  applied  to  the 
sample.  The  ."uclear  magnetization  is  tipped  by  180  degrees,  from 
its  equilibrium  value  .Mo  to  -Mo.  Following  the  pulse  the  nuclear 
magnetization  begins  to  return  toward  equilibrium.  After  a  time 
tau  (the  pulse  delay  time)  has  elapsed,  a  90-degree  pulse  is 
applied.  The  magnitude  of  the  FID  following  this  pulse  is  pro¬ 
portional  to  the  state  of  the  magnetization  at  time  tau.  The 
initial  magnitude  of  the  FID  is  recorded,  and  the  experime.nt  is 
repeated  for  a  number  of  different  pulse  delay  ti.mes.  A  plot  of 
the  FID  magnitude  versus  tau  provides  t.he  characteristic  time  for 
the  return  of  the  nuclear  magnetization  to  equilibrium.  This 
characteristic  time  is  Tl ,  the  spin-lattice  relaxation  time. 

Relaxation  times  in  NMR  are  influenced  strongly  by  molecular 
motions  which  alter  the  angles  and  distances  between  nuclei.  In 
many  compounds,  the  rate  and  nature  of  molecular  motion  com¬ 
pletely  determines  the  spin-lattice  relaxation  time  Tl .  These 
effects  are  discussed  in  Chapter  III,  following  a  discussion  of 
the  spin-spin  relaxation  time  and  the  NMR  solid  echo. 


F.  Spin-Spin  Relaxation 

One  effect  of  the  rf  pulse  is  to  align  the  nuclear  magnets  and, 
at  the  end  of  a  90-degree  pulse,  all  the  nuclei  are  precessing  in 
phase.  They  do  not  process  in  phase  forever,  though.  Because 
neignbori.ng  nuclei  are  magnetic,  each  nucleus  is'  found  in  a 
magnetic  field  that  is  the  sum  of  the  externally-applied  field 
and  a  local  field  due  to  the  presence  of  neighboring  nuclei. 
Hence,  there  is  a  spectrum  of  precession  frequencies  for  nuclei 
within  the  sample. 

If  the  spectrum  of  precession  frequencies  is  narrow,  then  all  t.ne 
nuclei  process  at  essentially  t.he  same  frequency.  Following  an 
•  rf  pulse,  they  remain  in  phase  for  a  long  ti.me,  i.e.,  m>any 

precession  periods.  If  the  spectrum  is  broad,  t.nen  there  are 
large  differences  in  tne  precession  rates  between  nuclei.  Fol¬ 
lowing  an  rf  pulse,  tne  phase  coherence  of  t.he  nuclei  is  quickly 
destroyed  by  tnese  differences  i.n  precession  frequencies. 

T.he  character isti c  ti.me  that  the  nuclei  precess  in  pnase,  or  t.he 
characteristic  ti.me  for  the  syste.m  of  nuclei  to  come  to  internal 
equi  1  ibr  i  u.m ,  is  called  t.ne  spin-spin  relaxation  ti.me,  T2  .  Since 
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Figure 
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the  nuclei  must  be  at  internal  equilibrium  before  they  can  be  at 
complete  equilibrium,  then  it  is  obvious  that  Tl ,  for  a  given 
material,  is  always  greater  than  or  equal  to  T2.  For  protons  in 
common  substances  in  typical  magnetic  fields,  T2  ranges  from  tens 
of  microseconds  to  hundreds  of  seconds. 

In  many  materials,  and  especially  in  solids,  experimental 
estimates  of  T2  are  relatively  easy  to  make.  First,  the  maximum 
magnitude  of  the  FID  is  noted.  The  FID  is  obscured  near  the  end 
of  the  rf  pulse  due  to  the  finite  recovery  time  of  the  receiver, 
so  the  FID  magnitude  is  recorded  just  after  the  receiver  system 
recovers.  Then  the  FID  response  is  followed  until  its  magnitude 
falls  to  (1/e)  of  its  former  value.  The  time  for  this  decay  is 
approximately  equal  to  T2 ,  though,  in  fact,  the  true  value  of  T2 
depe.nds  upon  the  exact  shape  of  the  FID. 

When  NMR  samples  contain  two  or  more  substances,  it  is  comm.on  to 
observe  more  than  one  T2  time.  A  typical  example  is  a  solid 
compound  that  has  absorbed  water;  one  T2  time  is  observed  for  t.he 
material  itself,  and  one  for  the  absorbed  water. 

The  relaxation  times  Tl  and  T2  are  important  because  many  explo¬ 
sive  compounds  possess  almost-unique  values  of  t.hese  parameters. 
Therefore,  an  NMR  'system  may  be  designed  to  identify  these 
compounds  in  the  presence  of  orher,  benign  materials. 


G,  Solid  Echo 

Most  solid  substances  exhibit  very  short  T2  times  and,  as  a  con¬ 
sequence,  much  of  the  FID  signal  is  obscured  by  the  dead  time  of 
the  NMR  system  (which  is  sometimes  called  the  recovery  time).  A 
two  pulse  technique,  producing  a  response  called  a  solid  echo, 
can  be  used  in  many  cases  to  recover  some  of  this  lost  signal. 

The  solid  echo  sequence  is  shown  in  Figure  II-4.  At  time  zero,  a 
90-degree  pulse  is  applied  to  the  nuclear  spin  system.  Following 
this  pulse,  the  spin  system  responds  with  an  FID,  but  much  of  the 
FID  signal  is  obscured  by  the  receiver  dead  time. 

Shortly  after  the  receiver  recovers,  as  shown  at  time  tau  in 
Figure  II-4,  a  second  90-degree  pulse  is  applied.  The  rf  phase 
of  this  pulse  is  shifted  by  90  degrees  from  that  of  the  first 
pulse.  The  effect  of  the  second  pulse  is  to  produce  an  NMR 
response,  called  a  solid  echo,  which  builds  in  time,  goes  through 
a  maximum,  and  then  decays  to  zero.  A  portion  of  the  solid  echo 
is  also  obscured  by  the  dead  time  of  t.he  receiver;  but  if  rhe 
pulse  spacings  are  properly  chosen,  then  the  maxi.mum  in  the  solid 
echo  will  be  observed,  as  shown  at  time  2 (tau)  in  Figure  I 1-4. 

The  solid  echo  can  be  used  to  improve  the  signal -to-noise  ratio 
of  t.he  system  by  improving  the  signal  response,  and  it  can  also 
be  used  to  determine  the  complete  shape  of  the  NMR  free  induction 
decay.  It  is  anticipated  that  the  solid  echo  technique  will  be 
used  in  any  single-sided  NMR  system  for  remote  detection. 


(dead  times) 


Solid  Echo 
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H.  Signal-to-Noise  Ratio 

The  fundamental  limit  of  sensitivity  for  an  NMR  system  is  deter¬ 
mined  by  the  signal -to-noise  ratio  (SNR)  of  that  system.  In 
evaluating  the  technology,  it  is  important  to  be  able  to  predict 
the  changes  in  SNR  that  will  result  from  changes  in  other  para¬ 
meters  . 


Laboratory  Systems 

For  proton  NMR  systems  in  typical  laboratory  configurations,  over 
a  wide  range  of  operating  parameters,  the  signal -to-noise  ratio 
scales  according  to  the  following  expression  (12), 


SNR  •(  f  (  Wo^QV/B) 


(II  -  1) 


where  f  is  the  filling  factor  of  the  sample  circuit  (i.e.,  the 
ratio  of  sample  to  sample  coil  volumes).  Wo  is  the  Larmor 
frequency,  Q  is  the  quality  factor  of  the  sample  circuit,  V  is 
the  volume  of  the  sample  coil,  and  B  is  the  bandwidth  of  the 
receiver/detector  circuitry. 

Equation  II-l  predicts  that  the  SNR  increases  as  the  Larmor 
frequency  to  the  (3/2)  power.  In  fact,  this  is  a  conservative 
prediction.  A  comprehensive  analysis  of  the  problem  (13)  re¬ 
sulted  in  the  conclusion  that  the  signal -to-noise  ratio  varies  as 
the  Larmor  frequency  to  the  (7/4)  power.  However,  over  the  range 
of  frequencies  under  consideration  here,  the  difference  between  a 
(3/2)  and  a  (7/4)  power  dependence  is  not  significant. 


Remote  Detection  Systems 

In  remote  detection  NMR,  the  precessing  nuclear  magnetization 
induces  a  signal  voltage,  i.e.,  the  FID,  into  a  detector  coil. 
The  voltage  induced  in  the  detector  coil,  Es ,  is  given  by  (6) 


Ss  =  NQAWoMo/Z^ 


(II  -  2) 


In  this  expression,  N  is  the  number  of  turns  in  the  detector  coil 
with  quality  factor  Q  and  cross-sectional  area  A.  In  addition. 

Mo  is  the  total  nuclear  magnetization  precessing  at  the  Larmor 
frequency  Wo,  and  Z  is  t.he  center-to-center  distance  between  the 
remote  region  of  precessing  magnetization  and  the  detector  coil. 
It  is  assumed  that  the  detector  coil  is  aligned  in  order  to 
maximize  the  induced  voltage  Es . 
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The  noise  voltage  for  a  remote  detection  system  can  be  estimated 
from  the  well-known  expression  for  Johnson  noise.  The  noise 
voltage  present  in  the  sample  circuit  is  given  by 


En  =  (  4kTBR) 


(II  -  3) 


where  R  is  the  resistive  component  of  the  impedance  of  the  sample 
circuit,  T  is  the  temperature  of  the  sample  circuit,  B  is  the 
bandwidth  of  the  system,  and  k  is  Boltzmann's  constant. 

Combining  Equations  II-2  and  II-3  provides  an  expression  useful 
for  estimating  the  SNR  of  prospective  remote  NMR  syste.ms .  The 
voltage  signal-to-noise  ratio  may  be  written 


SNR  =  (Es/En)  *  F 


(II  -  4) 


where  F  is  a  term  that  incorporates  all  of  the  factors  which  may 
affect  the  SNR  that  are  not  included  in  Equations  II-2  and  II-3. 
Examples  of  these  factors  include  the  effects  of  the  system 
recovery  time  and  the  FID  decay  time  T2 ,  the  effect  of  the  length 
of  t.he  rf  pulse,  the  effect  of  the  homogeneity  of  the  rf  field 
over  the  remote  region,  the  effects  of  the  precise  geometry  and 
orientation  of  the  sample  coil,  and  the  effects  of  noise  intro¬ 
duced  by  subsequent  amplification  of  the  NMR  signal.  In  general, 
F  is  difficult  to  determine  quantitatively.  However,  under  a 
fairly  wide  range  of  conditions,  qualitative  estimates  of  F  are 
relatively  easy  to  .make. 

Recalling  that  .Mo  =  Mv  *  Vm,  where  Vm  is  the  volume  of  the  remote 
region  and  Mv  is  the  nuclear  magnetization  per  unit  volume,  gives 


SNR  =  ( lo"^ ) ( 5/2k ) *  (WoMvVm/Z^)  * 


{ NAQ ) ( TBR 


-1/2  * 


(II  -  5 


T.he  first  term  in  Equation  II-5  is  a  constant;  the  10  ^  factor 
IS  included  to  allow' use  of  the  convenient  units  discussed  below. 
The  second  term  is  determined  by  the  particular  physical 
configuration  of  the  remote  NMR  syste.m.  That  is,  a  given  remote 
NMR  system  will  result  in  a  given  set  of  t.ne  parameters  Wo,  Vm, 
and  R;  and  Mv  will  be  determined  by  the  characteristics  of  t.he 
material  in  the  remote  region. 
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The  third  term  in  Equation  II-5  is  determined  primarily  by  the 
detector  coil  and  its  associated  circuitry.  If  the  sample  cir¬ 
cuit  is  parallel  tuned,  rather  than  series  tuned,  the  resistive 
component  of  the  circuit  impedance  is  given  by  R  =  QX,  where  X  is 
the  inductive  reactance  of  the  detector  coil.  With  this  substi¬ 
tution,  the  third  term  becomes 


Term  3 


(NAQ)  (TBR)'^'^^  =  (NA)  (Q/TBX) 


Utilizing  X  =  WoL,  where  Wo  is  the  Larmor  frequency  and  L  is  the 
inductance  of  the  sample  coil,  and  the  fact  that  the  maximum  Q  of 
the  sample  circuit,  Qm,  is  limited  by  the  system  bandwidth  to  the 
value  Qm  =  Wo/2TrB,  gives 


Term  3  =  (  NA/B  )  ( 1/2-frTL) 


Over  a  fairly  wide  range  of  parameters,  the  inductance  of  the 
detector  coil  is  expressed  by  (14) 


L  =  CN^d 


where  L  is  given  in  Henries,  N  is  the  number  of  turns,  and  d  is 
the  diameter  of  the  coil  in  centimeters.  The  constant  C  is  given 
by 


C  =  ( 10"^/2.54)/( 18  +  40(^/d)) 


(II  -  <^) 


where  i  is  the  length  of  the  coil  in  centimeters.  This  expres¬ 
sion  for  C  is  valid  for  iS/d)  ratios  greater  than  0.4,  and  it  is 
unlikely  in  remote  detection  NMR  that  an  (//d)  ratio  greater  than 
about  2.0  would  be  encountered.  Therefore,  the  constant  C  ranges 
from 


4.017  X  10 


-9 


<  C  <  1.158  X  10 


-8 
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Incorporating  the  fact  that  the  area  of  the  coil  A  =  TTd  /4, 
gives 


Term  3  =  (  1/4B  )  (rd^/2TC) 


Substituting  this  expression  into  Equation  II-5,  and  collecting 
the  numerical  constants'  into  the  first  term,  gives  the  following 
equation  for  the  signal -to-noise  ratio. 


SNR  =  (  10~^/8  )  (  5ir/k)  *  (WoMvVm/Z^)  * 


1/3)  (d-'/TC) 


3  * 


(II  -  7 


The  numerical  values  and  units  for  the  terms  used  in  Equation 
II-7  are  given  below. 


k  =  Boltzmann's  Constant, 

=  1.38  X  lO”^^  ergs/Kelvin; 

Wo  =  Larmor  frequency,  in  radians/second, 

=  2.675  X  10^  at  Ho  =  1,000  Gauss; 

Mv  =  Nuclear  magnetization  per  unit  volume 
in  the  sensitive  region,  in  abamp/cm, 

=  3.170  X  10  (for  6.6  x  10^^  spins/cc. 

Ho  =  1,000  Gauss,  and  Sample  Temp  =  300  K); 

Vm  =  Volume  of  sensitive  region,  in  cm^ ; 

Z  =  Center-to-center  distance  between  pickup 
coil  and  sensitive  region,  in  cm; 

d  =  Diameter  of  pickup  coil,  in  cm; 


B 


Bandwidth  of  the  detection  system,  in  Hertz; 


1 


T  =  The  temperature  of  the  sample  circuit  source 
resistance,  in  degrees  Kelvin; 


C  =  A  unitless  constant  that  depends  upon  the 
geometry  of  the  pickup  coil,  given  by 
Equation  II-6; 


F  =  A  unitless  term  incorporating  all  of  the 
factors  that  affect  the  SNR  not  otherwise 
included  in  Equation  II-7. 
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III.  NMR  AND  THE  MATERIALS  OF  INTEREST 


A.  Introduction 

Differences  in  the  NMR  relaxation  properties  of  materials  form 
the  basis  for  NMR  methods  of  explosives  detection.  Therefore,  a 
general  appreciation  of  the  reasons  why  different  materials 
exhibit  different  NMR  properties  is  central  to  an  understanding 
of  the  detection  technique.  In  the  following  discussion,  all 
comments  and  quantities  apply  specifically  to  proton  NMR,  and  not 
to  NMR  characteristics  of  other  nuclei. 

In  general,  the  relaxation  time  T1  for  protons  is  influenced  by: 
(a)  the  strength  of  the  magnetic  interaction  between  a  given 
proton  and  its  neighbors;  and,  (b)  the  relative  motion  between 
the  proton  and  its  neighbors.  When  the  frequency  of  this 
relative  motion  is  at,  or  very  near,  the  Larmor  frequency,  then 
T1  is  short.  When  the  frequency  of  this  motion  is  either  muc.h 
higher  or  much  lower  than  the  Larmor  frequency,  then  T1  is  long. 
The  strength  of  the  magnetic  interaction  depends  primarily  upon 
the  nearest  neighbor  distances  between  nuclei  and  upon  the 
presence,  or  absence,  of  paramagnetic  impurities. 

If  the  molecular  motion  near  the  Larmor  frequency  is  simple, 
i.e.,  if  there  is  only  a  single  type  of  motion  and  all  nuclei 
participate,  then  the  behavior  of  T1  is  relatively  easy  to 
characterize.  The  behavior  of  T1  is  generally  discussed  in  terms 
of  t,  the  correlation  time  for  the  molecular  motion.  If,  for 
example,  the  molecular  motion  is  translational  diffusion,  then 
the  correlation  time,  t,  can  be  thought  of  as  the  time  between 
diffusional  jumps.  Or,  if  the  motion  is  molecular  reorientation, 
then  the  correlation  time  can  be  thought  of  as  the  average  time 
it  takes  for  a  molecule  to  rotate  through  an  angle  of  one  radian. 
In  a  broad  sense,  the  correlation  time  may  be  thought  of  as  the 
inverse  of  the  frequency  of  the  molecular  motion. 

Figure  III-l  is  a  log-log  plot  of  T1  versus  the  correlation  time 
t  for  three  different  Larmor  frequencies  (A,  B  and  C).  In  the 
region  of  short  correlation  times,  i.e.,  for  high  temperatures  in 
most  substances,  the  frequency  of  molecular  motion  is  much  higher 
than  the  Larmor  frequency.  Here,  the  Tl  of  the  material  is  long 
and  increases  as  the  correlation  time  decreases.  Also,  in  this 
region,  Tl  is  independent  of  the  Larmor  frequency  chosen. 

When  the  frequency  of  molecular  motion  is  roughly  equal  to  the 
Larmor  frequency,  Tl  goes  through  a  minimum.  Since  Larmor  fre¬ 
quencies  can  be  altered  experimentally  by  changing  Ho,  separate 
Tl  minima  may  be  observed  and  these  are  shown  in  Figure  III-l. 

When  the  correlation  times  are  long,  i.e.,  at  low  temperatures, 
where  the  frequency  of  molecular  motion  is  low  compared  to  the 
Larmor  frequency,  then  Tl  is  again  long  and  lengthens  with 
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Molecular  Motion 
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increasing  correlation  time.  However,  in  contrast  to  the  short 
correlation  time  region,  the  specific  value  of  Tl  in  this  region 
depends  upon  the  specific  Larmor  frequency  chosen  for  observa¬ 
tion.  Therefore,  some  indication  of  t.he  nature  of  the  molecular 
motion  can  be  obtained  through  studies  of  the  Larmor  frequency 
dependence  of  Tl  relaxation  times. 

The  relaxation  time  T2  is  a  rather  complex  function  of  the  two 
parameters  that  influence  Tl,  i.e.,  the  frequency  of  molecular 
motion  and  the  strength  of  the  internuclear  coupling.  However, 
in  general,  T2  is  long  if  the  coupling  between  nuclei  is  weak, 
and  short  if  the  coupling  is  strong.  Both  Tl  and  T2  times  are 
influenced  by  changes  in  the  temperature  of  the  substance,  since 
temperature  changes  affect  the  frequency  of  molecular  motion.  Tl 
and  T2  are  also  sensitive  to  the  choice  of  Larmor  frequency, 
since  a  Larmor  frequency  near  a  frequency  of  molecular  motion 
results  in  a  short  Tl. 

Of  all  the  possible  molecular  motions  occurrina in  common  sub¬ 
stances,  three  often  occur  at  rates  close  enough  to  t.he  Larm.or 
freque.ncy  to  influence  Tl  and  T2 .  These  three  motions  are: 

(a)  molecular  tumbling,  (b)  molecular  diffusion,  and  (c)  re- 
orie.htation  of  mobile  groups  internal  to  the  molecule.  One,  or 
more,  of  these  three  motions  generally  determines  the  observed 
values  of  the  relaxation  times. 

In  liquids,  molecular  diffusion,  tumbling,  and  possible  internal 
motions  usually  proceed  at  rates  fast  compared  to  the  Larmor 
frequency  and,  therefore,  Tl  is  usually  long  (i.e.,  longer  than, 
say,  500  milliseconds).  In  addition,  since  the  molecules  are  in 
rapid  motion  relative  to  one  another,  the  coupling  between  nuclei 
is  weak  and  T2  is  also  long.  In  fact,  in  many  simple  liquids,  Tl 
and  T2  are  often  equal. 

In  solids,  molecular  diffusion  and  tumbling  usually  proceed  at 
rates  slow  compared  to  the  Larmor  frequency.  Therefore,  in  rigid 
crystalline  solids,  Tl  is  generally  long.  In  addition,  since  the 
atoms  are  fixed  in  place  and  the  nuclei  are  essentially  station¬ 
ary,  the  coupling  between  nuclei  is  strong  and  T2  is  short  (15). 

However,  most  solids  are  not  rigid  crystalline  materials.  Most 
common  substances  are,  in  fact,  composed  of  complex  molecules 
that  possess  side  chains  capable  of  a  wide  variety  of  internal 
motions.  If  many  side  chains  are  present,  then  it  is  likely  t.hat 
some  internal  motions  will  take  place  near,  or  above,  the  Larmor 
frequency  (16).  And,  of  course,  the  longer  the  side  chains,  the 
more  likely  it  is  t.hat  some  internal  motions  will  occur  near  the 
Larmor  frequency.  Therefore,  Tl  in  most  common  solid  substances 
IS  snort  (I.e.,  less  than  about  500  milliseconds)  and,  since  T2 
IS  always  less  than  or  equal  to  Tl ,  T2  is  also  short. 
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B.  Mixtures  of  Compounds 

In  use,  an  NMR  detection  system  will  often  encounter  many 
different  materials  simultaneously.  Therefore,  it  is  important 
to  determine  the  response  of  the  system  in  the  presence  of  two  or 
more  compounds. 

The  NMR  relaxation  time,  or  times,  exhibited  by  a  mixture  of 
compounds  may  be  difficult  to  analyze,  or  to  predict,  quanti¬ 
tatively.  The  results  depend,  in  large  part,  upon  whether  or  not 
the  molecules  of  the  compounds  are  in  intimate  contact.  If  not, 
and  even  if  the  compounds  are  in  very  close  proximity,  t.hen  the 
situation  is  described  by  the  separate  responses  from  each  of  the 
components  of  the  mixture.  If  so,  then  the  situation  is 
described  by  a  somewhat  more  complex  combined  response.  Each  of 
these  two  situations  is  discussed  below. 


Separate  Compounds 

The  description  of  the  relaxation  process  when  a  number  of 
separate,  noninteracting  compounds  are  involved  is  fairly 
straightforward.  The  case  of  two  separate  compounds  is  outlined 
here,  and  these  results  are  readily  generalized  to  account  for 
additional  compounds. 

Assume  that  the  two  compounds,  A  and  B,  are  allowed  to  come  to 
equilibrium  in  the  static  magnetic  field  and  a  180  degree  pulse 
is  applied.  The  nuclear  magnetization,  which  was  Mo  before  the 
pulse,  is  rotated  to  -Mo  immediately  afterward  (see  Chapter  II, 
Section  E).  At  time  t  later,  a  90  degree  pulse  is  applied  and 
the  Z-component  of  the  magnetization,  Mz ( t ) ,  is  determined  from 
the  FID  height. 

The  return  to  equilibrium  of  the  Z-component  of  the  magnetization 
is  given  by  the  following  expression: 


{.Mo  -  Mz{t))/2Mo  =  +  Be"^'^'^^  ,  (III  -  1) 


where  A  is  the  mole  fraction  of  compound  of  compound  A,  B  is  the 
mole  fraction  of  compound  B,  T1(A)  is  the  spin-lattice  relaxation 
time  of  compound  A,  and  T1{B)  is  the  spin-lattice  relaxation  time 
of  compound  B. 

Although  it  is  difficult  to  solve  Equation  III-l  in  general,  the 
solution  IS  easy  to  determine  i'n  specific  regions.  For  example, 
when  A  is  much  greater  than  B,  then  the  observed  T1  of  the 
separate  compounds,  T1{SC),  is  approximately  equal  to  Tl(A). 
Similarly,  when  B  is  much  greater  than  A,  then  Tl{SC)  is  appro.xi- 
mately  equal  to  T1(B). 
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It  IS  important  to  note,  however,  that  an  NMR  explosive  detection 
instrument  based  upon  Tl  differences  will  not  respond  to  the 
(l,e)  tim.e,  or  Tl  time,  calculated  from  Equation  Ill-i  for  a 
given  mixture.  Rather,  the  instrument  will  respond  to  the 
separate  FID  magnitudes  and  Tl  times  of  the  individual  com¬ 
ponents  of  the  mixture  falling  wit.hin  the  instrument  response 
window . 


Mixed  Compounds 

If  t.he  molecules  of  the  two  compounds  are  in  intimate  contact,  as 
in  solid  solutions,  then  t.he  state  of  the  relaxation  process  in 
one  compound  is  transmitted  to  the  ot.ner,  and  vice  versa.  This 
communication  takes  place  in  a  time  on  the  order  of  T2 .  In  the 
limit  where  T2  is  short,  i.e.,  the  communication  is  fast,  tnen 
the  resulting  relaxation  rate  is  given  by: 


1  TKMC)  =  A/TKA)  +  B/TKB), 


^  III  -  2) 


wnere  Ti(MC)  is  the  relaxation  time  for  the  mixture  and  the  other 
symbols  have  the  same  meanings  as  in  Equation  III-l. 

In  order  for  Equation  III-2  to  apply,  the  two  compounds  must  be 
mixed  on  a  molecular  level.  Therefore,  it  is  extremely  unlikely 
that  Tl(.MC),  from  Equation  III-2,  would  accurately  describe  the 
relaxation  observed  in  the  simple  mixture  of,  for  example,  two 
powders , 

It  is  not  out  of  the  question,  however,  that  situations  could 
arise  which  would  require  a  combination  of  Equations  III-l  and 
III-2  for  an  accurate  description.  Quantitative  analysis  of 
these  situations  would  be  difficult. 


C.  Explosive  Compounds 


Many  explosive  compounds  are  structurally  rigid  and  as  a  result, 
exhibit  long  proton  relaxation  times.  In  addition,  many  explo¬ 
sives  contain  nitrogen,  and  the  nitrogen  atoms  are  often  situated 
in  close  proximity  to  hydrogen  atoms.  As  a  consequence,  e.xplo- 
sive  compounds  often  exhibit  a  relaxation  effect  called  hydrogen 
nitrogen  level  crossing  whic.h  arises  from  an  interaction  between 
the  two  types  of  nuclei. 

Nitrogen  nuclei  possess  a  quadrupole  moment  and  tend  to  align 
along  electric  field  gradients  in  the  material.  This  alignment 
produces  nuclear  quadrupole  resonance  ( NQR )  frequencies  which, 
under  certain  conditions,  may  be  equal  to  t.he  proton  resonance 
frequency.  By  adjusting  the  Larmor  frequency,  or  magnetic  field, 
the  hydrogen  and  nitrogen  nuclei  can  be  brought  into  in  tnermal 
contact  and  can  exchange  energy.  This  is  called  level  crossi.tg. 

The  Tl  of  the  nitrogen  nuclei  is  generally  much  shorter  than  the 
Tl  of  the  hydrogen  nuclei.  The  energy  exchange  at  t.he  level 
crossi.ng  frequency  shortens  the  proton  Tl.  The  result  is  that 
many  explosive  compounds  exhibit  short  proton  Tl  times  at  certain 
specific  Larmor  frequencies  and  long  Tl  times  elsewhere. 

The  che.mical  properties  and  structures  of  a  variety  of  explosive 
compounds  are  discussed  below.  If  NMR  relaxation  times  for  these 
compounds  were  available,  they  are  included  along  with  comments 
on  i.mportant  considerations  related  to  the  NMR  detection  process. 

Black  Powder 

Black  powder  is  a  simple  mixture  of  potassium  nitrate,  charcoal 
and  sulfer.  It  contains  little  or  no  hydrogen.  Therefore,  black 
powder  is  not  detectable  by  proton  NMR. 


Compound  B 

Compound  B  is  a  cast  explosive  (18).  That  is,  it  is  cast  from  a 
mixture  into  predetermined  shapes.  Compound  B  is  an  inti.mate 
mixture  of  approximately  40%  TNT  and  60%  RDX.  Some  formulations 
may  include  up  to  63%  RDX  and  1%  of  an  unspecified  wax. 

It  Should  be  possible  to  estimate  tne  NMR  properties  zt  Compound 
3  from  t.he  NMR  properties  of  its  major  constituents.  The  k.nown 
NMR  properties  of  TNT  and  RDX  are  give.n  below,  and  tne  guidelines 
for  estimating  the  relaxation  properties  of  mixtures  are  iis- 
cussed  in  Section  B  of  this  chapter. 


Compound  C-4 


Compound  C-4,  sometimes  called  Harrisite  and  more  often  oy  tne 
generic  name  "plastic  explosive",  is  a  simple  mixture  of  RDX, 
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5.3%  Di ( 2-etny Ihexy 1 )  senacate,  2.1%  polyisobutylene,  and  1.6% 
motor  oil  (18).  The  NMR  properties  of  Di  (  2-ethylhe.xy  1  )  sebacate 
are  unknown.  The  NMR  properties  of  polyisobutylene  have  been 
recently  measured  (17)  and  this  rubbery  polymer  exhibits  a  rela¬ 
tively  short  T1  of  about  30  milliseconds  near  room  temperature. 
The  T1  of  motor  oil  is  unknown  but  is  likely  to  be  a  few  seconds 
or  less. 

Since  Compound  C-4  is  a  simple  mixture  of  these  four  compounds, 
Its  NMR  relaxation  properties  likely  approximate  those  of  RDX, 
Its  major  constituent. 


H.MX 


HMX,  or  1 , 3 , 5 , 7-tetranitro- 1 , 3 , 5 , 7-tetrazacyclooctane ,  is  a  pure 
explosive  compound  and  the  major  constituent  in  a  number  of 
plastic-bonded  explosives  (18).  The  chemical  structure  of  HMX  is 
shown  in  Figure  III-2. 

The  NMR  properties  of  HMX  are  not  available.  However,  the  chemi¬ 
cal  structure  suggests  that,  in  the  solid  form,  molecular  motions 
are  likely  to  be  slow.  Therefore,  Tl  is  expected  to  be  long.  In 
addition,  the  presence  of  a  large  number  of  nitrogen  atoms  indi¬ 
cates  that  significant  hydrogen/nitrogen  level  crossing  effects 
may  occur. 


PETN 


Pentaerythritol  tetranitrate ,  or  PETN,  is  a  white,  moldable 
compound  with  eight  hydrogen  and  four  nitrogen  atoms  per  molecule 
and  a  molecular  weight  of  316  (18).  The  chemical  structure  of 
PETN  is  shown  in  Figure  III-3. 

PETN  is  a  sensitive  and  very  powerful  explosive  with  a  detonation 
velocity  of  8,100  meters/second  (19).  It  is  used  extensively  in 
detonating  cords  and,  at  times,  as  a  primary  or  an  intermediary 
charge . 

PETN  exhibits  a  long  Tl  at  Larmor  frequencies  considered  suitable 
for  NMR  detection  (2).  At  3  MHz,  the  proton  Tl  is  approximately 
60  seconds.  In  common  with- many  explosives,  PETN  also  exhibits 
hydrogen/nitrogen  level  crossing.  Level  crossing  occurs  at  pro¬ 
ton  Larmor  frequencies  of  0.50  and  0.90  MHz  which  correspond  to 
applied  magnetic  fields  of  117  and  211  Gauss,  respectively. 


RDX 

RDX,  or  I , 3 , 5-trinitro-l , 3 , 5-triazocyclohexane ,  is  a  whit? 
moldable  compound  wit.h  six  hydrogen  and  six  nitrogen  atoms  per 
molecule  and  a  molecular  weight  of  222  (18).  RDX  is  the  major 
component  of  the  often-used  plastic  explosive  C-4.  The  chemical 
structure  of  RDX  is  shown  in  Figure  III-4. 
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1 ,3.5.7-Tetranitro-i  .3,5,7-tetrazacyclooctane  (HMX) 

Figure  III-2 
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Pentaerythritol  tetranitrate 
(PETN) 


Figure  III-3 
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RDX,  originally  developed  as  a  medicine  (20),  exhibits  aoout  the 
same  explosive  power  as  PETN.  However,  RDX  is  considerably  more 
stable  chemical Iv  and  is  much  less  sensitive  to  mechanical  shock 
(  19  )  . 

RDX  exhibits  very  long  T1  times  at  Larmor  frequencies  considered 
suitable  for  XMR  detection.  3  MHz,  the  proton  T1  of  RDX  is 

312  seconds  at  3  MHz  (2). 

RDX  also  exhioits  a  hydrogen/ nitrogen  level  crossing  effect  which 
lowers  t.he  proton  T1  time  at  certain  values  of  the  applied  field. 
Level  crossing  occurs  (2)  at  fields  of  423,  798  and  1,221  Gauss, 
which  correspond  to  proton  Larmor  frequencies  of  1.8,  3.4  and  5.2 
-MHz,  respectively. 


TATB 


TATB,  or  1  ,  3  ,  5-triamino-2 , 4 , 6-trinitrobenze.he  ,  is  a  bright  yellow 
pure  explosive  compound  with  the  chemical  structure  indicated  in 
Figure  III-5  { 18 ) . 

The  proton  NMR  relaxation  properties  of  TATB  were  determined  by 
Garroway  and  Resing  (21).  At  a  Larmor  frequency  of  60  .MHz,  TATB 
exhibits  fairly  constant  proton  T1  times  of  about  30  seconds  over 
a  wide  range  of  temperatures  indicating  little  molecular  motion 
at  high  frequencies.  No  data  were  available  concerning  hydrogen/ 
nitrogen  level-crossing  effects  in  TATB. 


TNT 


TNT,  or  2,4,6-trinitrotoluene,  is  a  light  brown  compound  with 
five  hydrogen  and  three  nitrogen  atoms  per  molecule  and  a  molec¬ 
ular  weight  of  227  (18).  TNT  is  considerably  less  powerful,  and 
less  sensitive,  than  RDX  (19).  The  chemical  structure  of  TNT  is 
shown  in  Figure  III-6. 

Compared  to  PETN  and  RDX,  TNT  exhibits  shorter  T1  times  at  the 
Larmor  frequencies  appropriate  to  NMR  detection  (2).  At  3  MHz, 
t.he  proton  Tl  of  TNT  at  3  MHz  is  approximately  25  seconds. 

Level -crossing  fields  in  TNT  are  comparable  to  those  found  in 
PETN.  Hydrogen/nitrogen  level  crossing  occurs  at  applied  fields 
of  185  and  204  Gauss.  These  fields  correspond  to  proton  Larmor 
frequencies  of  0.79  and  0.87  MHz,  respectively. 
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D.  Benign  Substances 

The  molecular  structures  and  NMR  properties  of  some  of  the  benign 
substances  that  may  be  encountered  during  the  inspection  process 
are  discussed  below.  It  should  be  noted  that  few,  if  any,  of 
these  materials  will  exhibit  the  hydrogen/nitrogen  level  crossing 
effects  characteristic  of  many  explosives. 


Bakel ite 


Bakelite  is  a  thermosetting  phenol  formaldehyde  resin  often  used 
as  an  electrical  insulator.  The  chemical  structure  of  bakelite, 
which  appears  to  be  quite  rigid,  is  shown  in  Figure  III-7. 

Though  no  NMR  data  are  available,  the  bakelite  structure  suggests 
that  molecular  motions  are  restricted.  Therefore,  it  is  likely 
that  T1  times  in  bakelite  are  long. 


Glass 


Glass  IS  a  silicate  and  contai.ns  no  protons.  Therefore,  there  is 
no  proton  NMR  response  from  glass. 


Metals 


The  presence 
process.  At 
reduced.  At 


of  metals  will  interfere  with  the  NMR  inspection 
best,  the  sensitivity  of  the  NMR  instrument  may  be 
worst,  no  signal  at  all  will  be  received. 


For  example,  if  an  explosive  is  placed  in  a  metal  container  or 
wrapped  in  metal  foil,  there  will  be  no  response  from  the  NMR 
instrument.  On  the  other  hand,  the  presence  of  metal  items,  or 
of  metal  foil  that  does  not  enclose  the  explosive,  will  produce  a 
loss  of  sensitivity  that  m.ay  or  may  not  be  nominal  ,  depending 
upon  the  relative  placement  and  orientation  of  the  materials. 


Plexiglas 


Polymethyl  methylacrylate  is  a  hard  clear  thermoplastic  polymer 
with  trade  names  such  as  Plexiglas  and  Lucite.  The  structure  of 
polymethyl  methylacrylate  is  shown  in  Figure  III-8.  At  room  tem¬ 
perature  and  at  a  Larmor  frequency  of  19  MHz,  Plexiglas  exhibits 
a  Tl  of  88  milliseconds  (22). 

Polymethyl  methylacrylate  has  short  simple  side  chains  that 
terminate  in  methyl  groups.  This  is  a  molecular  structure 
capable  of  exhibiting  a  wide  range  of  internal  motion  frequencies 
in  the  solid  state,  and  these  motions  are,  no  doubt,  responsible 
for  tne  short  Tl  times  ooserved  for  this  compound. 


Bakelite 

(phenol-formaldehyde  resin) 


CHj 


Figure  III -7 
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Molecular  Structure  of 
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Molecular  Structure  of  Polyethylene 


Figure  III -8 
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Pc Ivetnv lene 


Polyethylene  is  a  soft  clear  thermoplastic  polymer  with  the 
structure  shown  in  Figure  III-8.  It  is  a  very  simple  linear 
polymer  with  a  repeating  group  consisting  of  one  carbon  and  two 
hydrogen  atoms.  A  single  Tl  measurement  at  room  temperature  and 
19  MHz  gave  92  milliseconds  (22). 

The  opportunities  for  complex  motions  are  limited  in  polyethyl¬ 
ene,  but  "crankshaft"  rotations  and  chain  end  diffusion  may  occur 
over  a  wide  range  of  frequencies.  The  presence  of  these  motions 
is  likely  the  reason  for  t.he  short  Tl  observed  in  this  compound. 


Rubber 


The  structure  of  pol y isoprene ,  or  natural  rubber,  is  shown  in 
Figure  III-9.  Tl  times  observed  in  rubber  are  quite  short,  e.g., 
34  milliseconds  at  19  MHz  at  room  temperature  (22). 


T.his  compound  possesses  a  simple  repeating  unit  with 
side  chains.  Both  the  side  chains  and  the  repeating 
capable  of  a  wide  range  of  motional  frequencies,  and 
consistent  with  the  experimental  observation  of  very 
times. 


very  simple 
units  are 
this  IS 
short  Tl 


Styrofoam 

The  principal,  and  perhaps  only,  component  of  styrofoam  is  poly¬ 
styrene.  The  structure  of  polystyrene  is  shown  in  Figure  III-IO. 
Tl  ti.mes  for  styrofoam  at  19  MHz  (22)  range  from  0.75  to  0.82 
seconds  which  are  somewhat  longer  t.han  those  observed  for  many 
other  polymer  compounds. 

Polystyrene  is  a  linear  polymer  with  a  carbon  backbone.  Phenyl 
groups  are  connected  to  every  other  carbon  on  the  chain.  The 
size  of  the  phenyl  groups  hinders  many  of  the  possible  modes  of 
.molecular  motion  and  likely  slows  down  side  chain  reorientations 
and  "crankshaft"  rotations  considerably. 


Vinyl  Materials 

Vi.nyl  polymers  are  used  to  fabricate  a  wide  variety  of  ite.ms.  A 
typical  example  is  electrical  tape,  which  is  an  adhesive-backed 
vinyl  compound.  The  measured  Tl  time  for  electrical  tape  is  122 
milliseconds  at  room  temperature  and  a  Larmor  frequency  of  19  MHz 
>'22;.  Compared  to  the  Tl  times  for  explosives,  t.he  relaxation 
times  for  vinyl  compounds  are  likely  to  be  very  short. 
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Wood  Paoer 


The  principal  constituent  of  wood  and  paper  is  cellulose. 
Hawortn  diagram  of  the  repeating  unit  of  a  cel lu lose  cnain 
shown  in  Figure  III-ll.  In  botn  wood  and  paper,  tnese  cha 
of  irregular  length  and  are  arranged  in  an  amorpnous  fasni 

Wood  and  paper  e.xhibit  snort  Tl  times  over  a  wide  ranae  cf 
frequencies  (22).  An  e.xamination  of  Figure  III-ll  indicat 
a  variety  of  molecular  motions  are  possible  in  cellulose, 
ing  side  chain  reorientation,  "crankshaft''  rotations  of  a 
of  repeati.ng  subunits,  and  e.nd-of -chai.n  diffusive  motions, 
likely  occur  over  a  wide  range  of  frequencies  and  give 
t.ne  observed  rela.xation  oehavior. 
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IV. 


SYSTEM  COMPONENTS  FOR  REMOTE  NMR 


A.  Introduction 

The  purpose  of  this  program  is  to  explore  and  evaluate  the 
potential  of  nuclear  magnetic  resonance  (NMR)  for  remote  detec¬ 
tion  applications.  Regions  that  are  accessible  from  only  one 
side  are  of  particular  interest. 

Si.ng le-s ided  NMR  remote  detection  systems  differ  from  laooratory 
NMR  systems  in  that  the  required  region  of  homogeneous  DC  field, 
i.e.,  Ho,  must  be  produced  remotely.  The  required  rf  field.  Hi, 
must  also  be  produced  remotely.  As  in  the  case  of  laboratory  NMR 
systems,  .Hi  in  remote  detection  systems  must  oe  directed  perpe.n- 
dicular  to  the  steady  field  Ho. 

The  following  section  of  this  chapter  presents  alternate  methods 
of  producing  remote  regions  of  DC  magnetic  field.  T.he  tnird 
section  of  this  chapter  discusses  alternate  methods  for  producing 
remote  regions  of  rf  field.  In  the  final  section,  a  figure-of- 
merit  is  defined  which  allows  comparison  of  alternate  systems  of 
NMR  remote  detection. 


B.  DC  Field  Configurations 

An  NMR  system  must  be  capable  of  producing  a  magnetic  field  that 
is  constant  in  time,  i.e.,  a  DC  magnetic  field.  In  most  cases, 
It  IS  desirable  that  this  DC  field  be  homogeneous  in  the  region 
to  be  examined.  In  the  case  of  a  single-sided  NMR  system,  this 
homogeneous  region  of  steady  magnetic  field  must  be  produced 
remotely. 

Several  alternate  methods  of  producing  remote  fields  are  dis¬ 
cussed  in  the  following  sections.  First,  the  toroidal  geometry 
of  Inside-Out  NMR  is  presented.  Then,  a  discussion  of  U-shaped 
iron  magnet  technology  is  followed  by  a  detailed  treatment  of 
opposed  coils  and  solenoids. 


Opposed  .Magnet  Configuration 

An  opposed  .magnet  configuration,  termed  Remote  Inside-Out  NMR, 
was  developed  as  a  tool  for  borehole  logging  '5,6,7'.  The' 
homogeneous  region  of  magnetic  field  is  produced  by  placing  two 
magnets  along  a  line,  as  shown  in  Figure  IV- 1  ‘7).  The  li.ne  is  a 
line  of  symmetry  for  the  syste.m,  and  t.he  poles  of  the  two  magnets 
are  opposed. 

T.his  configuration  results  in  a  toroidal  region  of  homogeneous 
field.  The  toroidal  region  of  homogeneous  field  is  produced  in 
t.he  plane  that  is  equidistant  between  the  two  magnets  and 


OPPOSED  MAGNETS  PRODUCE  REGION 
OF  HOMOGENEOUS  RADIAL  MAGNETIC 
FIELD  NEAR  MIDPLANE 


Figure  17-1 


perpendicular  to  the  line  connecting  tnem.  The  field  is  directed 
radially  outward,  perpendicular  to  the  line  of  symmetry. 

In  practice,  the  homogeneous  region  of  field  may  be  produced  cy 
permanent  magnets,  electromagnets,  or  superconducting  coils. 
However,  since  the  homogeneous  region  results  from  magnets  in 
opposition,  the  .maximum  field  strengt.h  is  limited.  This  is 
illustrated  by  the  data  (7)  shown  in  Figure  IV- 2. 

In  t.his  example,  an  electromagnet  capable  of  a  maxi.mum  field 
strength  of  about  15  Kilogauss  in  a  typical  laooratory 
configuration  produced  a  .homogeneous  field  of  only  about  120 
Gauss  in  the  opposed  magnet  configuration.  In  the  laboratory 
configuration,  the  field  was  produced  over  a  region  10  cm  in 
diameter  and  2.5  cm  long,  i.e.,  a  2.5  cm  pole  face  gap.  In  the 
opposed  magnet  configuration,  the  field  was  produced  over  a 
toroid  approximately  14.5  cm  in  diameter,  with  a  cross-sectional 
diameter  of  approxi.ma te  1  y  3  cm. 


T-Shaped  .Magnet  Conf igurat ion 

Magnets  with  U-shaped  yokes  can 
fields  in  remote  regions.  Both 
magnets  are  suitable.  A  sketc.h 
Figure  IV- 3. 


be  used  to  produce  magnetic 
electromagnets  and  permanent 
of  such  a  syste.m  is  shown  in 


The  U-shaped  yoke  configuration  has  been  recently  used  in  NMR 
systems  to  measure  moisture  levels  in  structural  concrete  and  in 
subsurface  soils  (3,4).  The  remote  region  of  magnetic  field  is 
not  as  well  defined  as  that  produced  by  the  opposed  magnet 
configuration  and,  in  addition,  is  not  as  homogeneous.  That  is, 
field  gradients  within,  and  near,  the  region  of  constant  field 
are  greater  in  the  U-shaped  configuration  tnan  in  the  opposed 
magnet  configuration. 


I.n  fact,  t.he  cylindrical  region  of  uniform  field  sketched  in 
Figure  IV-3  represents  a  gross  simplification  of  the  actual 
situation.  Between  the  two  poles,  within  the  indicated  region  of 
uniform  field,  t.he  lines  of  constant  field  actually  conform  to 
tne  shape  of  a  saddle. 


The  .maximu.m  strengt.h  of  tne  remote  field  produced  oy  a  U-shaped 
magnet  configuration  would  naturally  be  expected  to  oe  somewhat 
greater  than  that  produced  by  configurations  of  opposed  magnets 
or  coils.  Available  information  on  the  system  used  to  measure 
.moisture  levels  in  concrete  indicates  that  this  is,  in  fact,  the 
case . 

The  inside  dimension  of  the  iron  yoke  for  the  magnet  used  in  tne 
system  for  concrete  .moisture  level  measurements  is  11  incnes 
2~.3  cm,.  At  maxi.mum  current,  t.his  system  produces  a  field  of 
494  Gauss  at  a  maximum  remote  dista.nce  of  6  incnes  :  15.2  cm ;  . 
This  distance  is  very  close  to  the  diameter  of  t.ne  remote  field 
toroid  in  the  opposed  magnet  example  discussed  aocve .  The  field 
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strength  in  the  remote  region,  though,  exceeds  that  of  tne 
opposed  magnet  example  by  more  than  a  factor  of  four. 


Opposed  Coil  and  Solenoid  Configurations 

Current-carrying  coils  and  solenoids  can  also  be  used  oo  produce 
remote  regions  of  homogeneous  .magnetic  field.  In  particular, 
configurations  of  coils  and  solenoids  designed  to  produce  fields 
that  partially  cancel  may  be  particularly  useful  12,23).  These 
configurations  are  commonly  called  opposed  coils  or  opposed 
solenoids . 

A  sc.he.matic  of  an  opposed  coil  pair  is  shown  in  Figure  IV- 4.  The 
system  consists  of  two  concentric  coils  of  different  radii.  In 
the  coordinate  system  used  here,  the  ce.nter  of  t.he  coils  is 
assumed  to  be  the  origin,  and  the  Z-axis  is  perpendicular  to  the 
plane  of  t.he  coil  windings.  The  currents  in  t.he  coils  flow  i.n 
opposition;  that  is,  they  produce  magnetic  fields  whicn  tend  to 
cancel  . 


Since  the  radii  of  the  two  coils  are  different,  t.he  fields  tnat 
they  produce  will  not  cancel  in  all  regions  of  space,  even  if 
tney  cancel  in  a  specific  region.  For  example,  assume  that  tne 
currents  are  adjusted  so  that  t.he  magnetic  field  at  t.he  origin  of 
tne  coordinate  system  is  zero.  Along  the  Z-axis,  the  field  due 
to  the  inner  coil  will  fall  off  more  rapidly  with  increasing  Z 
than  the  field  due  to  the  outer  coil.  Therefore,  as  Z  increases 
from  t.he  origin,  the  total  magnetic  field,  i.e.,  the  vector  sum 
of  t.he  fields  from  each  of  the  two  coils,  must  first  increase. 

However,  when  Z  becomes  large,  the  fields  from  eac.h  of  t.he  coils 
go  to  zero  and,  hence,  the  total  field  is  zero.  Therefore,  along 
the  positive  Z-axis,  the  magnetic  field  from  this  coil  pair  is 
zero  at  the  origin,  increases  with  increasing  Z,  goes  through  a 
.maxi.mum,  and  then  falls  back  to  zero. 

This  behavior  is  shown  in  the  plot  of  Figure  IV-5.  T.he  point  of 
maximum  field  occurs  at  a  distance  Zm  from  t.he  origin.  The 
precise  location  of  Zm  depends  upon  t.he  ratio  of  the  radii  of  the 
two  coils  and  t.he  ratio  of  the  two  current  flows.  In  this 
example,  t.he  radius  of  t.he  larger  coil  is  assumed  to  be  10  units 
and  t.hat  of  the  smaller  coil  5  units.  Each  produces  a  field  of 
unity  magnitude  at  the  origin.  The  maxi.mum  in  t.he  resultant 
field,  with  a  magnitude  of  0.369,  occurs  at  5.73  units  along  the 
Z-axis  . 


Aoout  the  poi.nt  Z.m,  the  homogeneity  of  the  .magnetic  field  can  be 
expressed  within  specified  limits.  T.he  size  of  tne  region  of 
ncmogeneous  field,  i.e.,  the  homogeneous  volume,  depends  not  only 
upon  t.ne  li.mits  specified,  out  also  upon  the  ratios  of  t.ne  coil 
radii  and  current  flows. 

T.he  t.nree  di.me.nsiona  1  shape  of  the  homogeneous  region  is  quite 
complex.  To  a  crude  approximation,  thougn ,  tne  homogeneous 
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region  can  be  assumed  to  be  spherical  in  shape.  This 
approximation  was  made  in  the  analyses  that  follow  in  a  later 
section . 

The  magnitude  of  the  field  at  its  maximum  and  the  size  of  the 
homogeneous  region  can  he  altered  by  offsetting  the  two  coils 
from  one  another.  This  is  shown  in  Figure  IV-6,  wnich  is  a 
sketch  of  opposed  coils  offset  by  a  distance  S.  The  origin, 
i.e.,  Z  =  zero,  is  assumed  to  be  at  the  center  of  the  coil 
closest  to  Zm,  the  point  of  the  remote  field  maximum.  The 
configuration  shown  in  Figure  IV-6  is  termed  a  negative  shift. 

If  S  extends  in  the  -Z  direction,  this  is  termed  a  positive 
shift . 

It  IS  of  interest  to  note  that  there  are  particular  combinations 
of  coil  radii,  coil  currents  and  offset  distances  that  produce  a 
region  around  the  point  Zm  where  not  only  the  first  derivative 
but  also  the  second  derivative  of  the  field  with  respect  to 
distance  goes  to  zero.  By  analogy  to  Helmholtz  coils,  these  are 
called  Inside-Out  Helmholtz  configurations  (23). 

Inside-Out  Helmholtz  coils  produce  a  larger  region  of  homogeneous 
field  than  do  pairs  of  coils  with  arbitrary  spacing  and  currents. 
A  sketch  of  an  Inside-Out  Helmholtz  pair  is  shown  in  Figure  IV-7, 
and  the  corresponding  plot  of  fields  along  the  Z-axis  is  shown  in 
Figure  IV-8.  Inside-Out  Helmholtz  opposed  coils  require  a 
positive  shift,  which  means  that  the  smallest  of  the  two  coils  is 
closest  to  the  remote  region  of  homogeneous  field. 

Fields  that  partially  cancel  can  also  be  generated  by  opposed 
solenoids.  As  shown  in  Figures  IV-9  and  IV-10,  opposed  solenoids 
can  be  either  aligned  or  offset.  These  configurations  are 
analagous  to  the  coil  configurations  discussed  above,  and  the 
coordinates  used  for  the  two  systems  are  identical. 


C.  RF  Coil  Configurations 

In  addition  to  the  steady  magnetic  field,  an  NMR  remote  detection 
system  must  also  produce  a  radio  frequency,  or  rf,  field  within 
the  region  of  interest.  This  rf  field  must  be  aligned  perpendi¬ 
cular  to  the  steady  field  in  order  to  produce  an  NMR  response. 

Because  of  practical  limitations  on  the  amount  of  rf  power  that 
can  be  conveniently  generated  and  delivered,  it  is  generally  not 
feasible  to  use  inefficient  rf  coil  configurations  (e.g.,  opposed 
loops)  when  short  rf  pulses  are  required.  A  consequence  of  this 
constraint  is  that  appropriate  rf  coil  configurations  will  likely 
produce  an  rf  field  that  is  inhomogeneous  over  the  useful  region 
of  DC  field.  However,  with  careful  design,  it  should  be  possible 
to  mi.nimize  the  effects  of  this  inhomogeneity. 

Three  rf  coil  configurations  are  examined  below.  A  spiral 
configuration  is  discussed  first,  and  this  is  followed  by 
discussions  of  the  solenoid  and  semi-toroid  configurations. 
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Schematic  of  Inside-Out  Helmholtz  Coils 
to  generate  the  DC  Magnetic  Field 


Figure  IV-7 


Figure  IV-10 


The  spiral,  or  spi  ra  1 -wound ,  coil  is  a  two  di.nens  icna  1  coil 
wnich,  for  analysis,  car.  be  considered  to  oe  a  series  of  con¬ 
centric  single-turn  coils.  Spiral  coils  are  someti.c.es  called 
spiral-wound  coils  or  pancake  coils.  A  spiral  coil  is  sncwn 
sche.Tiatica  1 1  y  in  Figure  17-11, a;. 

Near  tne  coil,  and  along  t.he  a.xis  of  symmetry,  a  spiral  coil 
produces  a  field  perpendicular  to  the  plane  of  the  coil  windinc. 
This  IS  illustrated  in  Figure  IV-ll(b). 

For  applications  as  an  rf  coil  in  NMR  re.mote  detection,  tne 
spiral  coil  must  be  used  in  con;^unction  wit.h  a  syste.m  t.nat 
produces  a  DC  field  perpendicular  to  the  symmetry  a.xis  of  t.ne 
spiral  coil.  Of  t.ne  syste.ms  discussed  aoove ,  only  the  U-shaped 
configuration  satisfies  t.nis  constrai.nt. 

There  is  one  obvious  drawback  to  the  use  of  a  spiral  rf  coil  in 
a  i.ng  le-s  ided  NMR,  which  is  t.hat  the  rf  field  is  not  focused,  or 
directed,  toward  tne  region  of  interest.  If  a  spiral  rf  coil  is 
placed  flush  agai.nst  t.ne  region  to  be  examined,  so  that  t.ne  rf 
field  penetrates  into  the  region,  t.hen  an  rf  field  of  equal  mag¬ 
nitude  IS  also  produced  outside  the  region  of  interest. 

This  inefficiency  is  reflected,  in  part,  by  the  rf  power 
requirements  of  the  system  designed  and  built  to  measure  moisture 
levels  in  structural  concrete  (3).  In  this  system,  an  rf  pulse 
power  of  200  kilowatts  is  required  to  produce  a  30  microsecond, 

90  degree  pulse  at  a  Larmor  frequency  of  2.1  .MHz. 

Solenoid  Coi 1 

A  solenoid  coil,  which  is  commonly  used  in  laboratory  NMR 
instruments  as  the  sample  coil,  is  shown  schematically  in  Figure 
17-12.  From  this  figure  it  is  obvious  that  the  solenoid  coil 
produces  a  remote,  or  fringe,  field  that  in  some  regions  is 
parallel  to  t.he  coil  axis  of  symmetry. 

The  direction  of  t.he  rf  field  produced  by  a  solenoid  coil  makes 
It  appropriate  for  use  wit.n  either  the  toroidal  DC  field  oroduced 
by  t.ne  opposed  magnet  geometry  or  the  localized  DC  field  produced 
oy  opposed  loops  such  as  t.ne  Inside-Out  Helmnoltz.  It  should  be 
pointed  out  t.hat  the  symmetry  of  t.ne  field  oroduced  bv  a  solenoid 
coil  IS  ideally  suited  to  the  toroidal  geometry  of  opposed  magnet 
NMR.  It  IS  less  t.nan  ideal  for  use  wit.n  a.n  I.nside-Cut  Hel.mholtz 
geometry  since  its  field  pattern  is  not  specifically  directed 
toward  t.ne  remote  region  of  homoce.neous  field. 


itioroid 


The  semitoroid  coil  '24)  can  be  hnoucnc  of  in  two  ways:  a  as  s 

solenoid  coil  which  has  been  oent  into  a  semicircle;  or,  o  as  a 
dcnut  snaped,  or  toroidal  coil,  which  has  oeen  cut  in  half.  A 
SKetcn  of  a  se.mitoroid  coil  and  its  acccmpanyi.ng  field  lines,  is 
snown  in  Figure  17-13. 


The  se.mitoroid  rf  coil  produces  an  Hi  field  which  is  perpendic¬ 
ular  to  the  axis  of  symmetry  of  the  coil,  i.e. ,  tne  Z  axis.  Th 
means  that  tna  semitoroid  coil  is  appropriate  for  use  witn  tne 
Inside-Cut  Helmholtz  configuration  for  generating  the  DC  field. 
In  addition,  tne  semitoroid  coil  produces  an  rf  field  pattern 
t.nat  15  directed  toward  the  remote  region  of  homogeneous  field. 
Therefore,  the  semitoroid  coil  should  be  more  efficie.nt  in 
single-sided  NMR  syste.ms  than  either  spiral  or  solenoid  coils. 
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D.  Figure-of -Merit 

In  order  to  evaluate  a.nd  compare  alternate  spectrometer  configu¬ 
rations,  a  figure-of -merit  is  required.  A  f igure-of-merit  allows 
potentially-interesting  geometries  to  be  screened  without  the 
necessity  of  performing  detailed  calculations  for  each  one. 

To  arrive  at  t.he  f  igure-of-merit ,  assume  that  an  NMK  spectro¬ 
meter  produces  a  remote  region  of  field  and,  further,  that  the 
center  of  this  remote  region  is  located  a  distance  Zm  from  the 
plane  of  the  coils  producing  the  field.  The  total  nuclear 
.magnetization,  .Mo,  wit.hin  this  remote  region  is  proportional  to 
tne  product  of  t.he  volume  of  the  remote  region,  Vm,  times  the  DC 
field  strengt.h.  Ho.  That  is. 


-Mo  M,  7m  *  Ho 


:7 


Tne  quantity  .Me  is  i  measure  of  t.ne  nuclear  magnetizati 
t..ted  oy  tne  spectrometer,  t.nat  is  availaole  to  oe  m.ani 
I  .n  sccie  se.nse,  .'.c  is  3  f  i  uure  — o  f  — m^er  1 1 .  Heweu'er  ,  .'■!o  is 
f -gure-c  f-.merit  needed  to  :.dentify  promising  spectrom'et 
fiuuraticns  for  re.mote  detecti'on. 
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Scnematic  of  5ern:-.oroic]  RF  Coil 

Figure  17-13 


£3  a<  Me  *  Wo  Z' 


_  tne  larmor  frequency .  ComDining  ill  and  'II',  asinc 
nae  Wo  is  propomcna’  zo  Ho,  and  ledoinc:  z  =  Z~  allows 
gore-of -men  t ,  to  oe  defined: 


?s  =  Vm  *  Ho".  Z;n''  '  IV  -  3  i 


wnere  ..o  i 
tne  fact  t 
F  3  f  -  n  0  z  1 


?s  is  a  meas'cre  of  tne  N’MR  signal  voltage  induced  in  an  rf  pickup 
coil  for  a  given  spect rorneter  configuration.  Si.nce  the  '.’olurnie  of 
tne  homogeneous  region  of  field  Vrn  is  proportional  to  Zm  cuced 
'at  least  to  a  first  approximation,  for  the  cases  studied  here', 
and  since  Ho  depends  only  upon  t.ne  DC  currents  applied  to  the 
coils,  It  follows  that  Fs  provides  a  sim.ple  .measure  for  comparing 
competing  configurations. 
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ri .  InrroGuct  icn 

Desirable  qualities  of  a  sincle-sided  NMR  system  for  remote 
detection  applications  include: 

*  a  penetration  depth  of  15  -  20  cm,  or  more; 

*  pcrtaoility,  perhaps  witn  a  mo'/eaole  inspection 
nead ;  and, 

*  a  well-defined  region  of  homogeneous  magnetic 
field  with  a  strength  of  at  least  1,000  Gauss. 

Tc  elacorate,  a  minimum  penetration  deptn  in  the  range  of  15  to 
2  0  c.m  Will  be  required  in  many  of  the  anticipated  inspection  ap¬ 
plications.  In  ge.neral,  of  course,  the  greater  the  penetration 
deptn,  tne  better. 

In  some  applications,  portability  of  t.he  resulting  system  may  be 
required;  in  others,  portability  may  be  very  desirable.  Trans¬ 
portability,  in  contrast  to  portability,  will  be  important  in 
virtually  all  potential  applications. 

•A  field  strength  of  at  least  1,000  Gauss  will  be  required  for  an 
adequate  signal -to-noise  ratio  and,  in  addition,  for  utilizing 
such  potentia 1 1 y-unique  NMR  signatures  as  hydrogen-  nitrogen  level 
crossing.  A  well-defined  region  of  homogeneous  field  will  ice 
i.mportant  in  applications  where  the  precise  location  of  the 
material  of  interest  must  be  determined. 

With  these  goals  in  mind,  five  alternate  syste.m  configurations 
for  remote  NMR  detection  were  considered.  These  five  system 
configurations  are  briefly  descrioed  below. 


1.  .q  syste.m  utilizing  the  toroidal  geometry  of  the  opposed 
magnet  co.nf  iguration .  I.n  t.nis  syste.m,  opposi.ng  magnets 
are  used  to  produce  t.ne  re.mote  region  of  nomoge.necus  DC 
field  and  a  solenoid  coil  is  used  to  produce  tne  rf  field. 

2.  An  instrum.ent  configuration  in  wnicn  the  DC 
produced  oy  a  T-snaped  iron-core  magnet  and 
-5  producec  by  a  flat  spiral  coil. 


- 


58 


4.  A  aysre.Ti  configuration  in  whicn  opposed  coils  or  sole¬ 
noids  are  used  to  produce  a  homogeneous  region  of  rf 
field,  and  the  required  DC  field  is  produced  oy  an  al¬ 
ternate  configuration,  e.g.,  a  semitoroid. 

5.  An  instrument  configuration  in  which  both  the  DC  and 

rf  fields  are  produced  by  semitoroids  and  are  relatively 
ncnuniform. 


Development  work  has  been  carried  out  previously  on  the  first  two 
configurations  and  the  results  of  this  work  are  described  in  tne 
next  section.  This  is  followed  by  a  discussion  of  advanced- 
design  single-sided  NMR  remote  detection  systems.  In  t.ne  final 
section  of  t.nis  chapter,  the  proposed  si.ngle-sided  systems  are 
evaluated  usi.ng  the  figure-of -merit  developed  in  C.hapter  IV. 


3.  Present  Technology 

.Pecent  work  in  re.mote  detection  NMR  nas  been  directed  toward  two 
pri.tcipal  applications.  Re.mote  Inside-Out  NMR  was  conceived  and 
developed  as  a  borehole  logging  tool  for  geopnysical  evaluation. 
Single-sided  NMR  using  a  U-shaped  magnet  was  developed  for 
-teasure.ments  of  moisture  in  structural  concrete  and  subsurface 
soils.  These  two  technologies  are  discussed  in  the  following 
sections . 


Inside-Out  NMR 

Inside-Out  Nf-lR  was  developed  for  t.he  re.mote  detection  of  fluids 
trapped  in  porous  rock  formations  adjacant  to  a  borenole  (5,6,7). 
An  Inside-Out  NMR  system  combines  two  of  the  compone.nts  dis¬ 
cussed  in  Section  IV. B:  an  opposed  magnet  system  (to  produce  a 
toroidal  region  of  DC  field)  and  a  solenoid  coil  (to  produce  a 
perpendicular  rf  field  in  the  remote  region). 

A  sketch  of  an  Inside-Out  NMR  system  is  shown  in  Figure  V- 1  (7). 
The  rf  coil  is  located  halfway  between,  and  on  the  axis  of,  the 
two  magnets  t.hat  produce  the  DC  field.  The  magnetic  field  in  the 
toroidal  region  is  directed  radially  outward,  and  the  fringe 
field  of  the  solenoid  rf  coil  i.ntersects  t.he  DC  field  at  a  90 
degree  a.ngle. 

When  an  rf  pulse  is  applied  to  the  sole.noid  coil,  as  shown  in 
Figure  V- 1  ,  t.he  nuclear  magnetization  in  the  re.mote  region  is 
tipped  from  its  equilibrium  position.  The  duration  of  the  rf 
pulse  and  the  strength  of  the  rf  field  determine  the  tip  angle. 

At  tne  end  of  the  rf  pulse,  the  .nuclear  magnetization  precesses 
about  tne  steady  field,  as  shown  i.n  Figure  V- 2  .  The 

precessing’ magnetization  induces  a  voltage  across  the  leads  of 
the  solenoid  coil,  i.e.,  tne  NMR  free  induction  decay  signal. 
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OSCILLATING  MAGNETIC  FIELD 
FROM  PRECESSING  NUCLEAR 
MAGNETIZATION 


While  tne  aeometry  or  Inside-Ouh  \'MR  is  ideal  for  rerr.coe  aeoec- 
tion  in  well  logging  applications,  it  is  less  suitaole  frr  incse 
applications  where  access  is  li.T.ited  to  only  one  side.  In 
addition,  at  the  penetration  deptns  required,  it  is  inli.<el  y  t.nat 
tne  Inside-Cut  geometry  could  produce  tne  required  field 
strencths  of  1,000  Gauss,  even  with  superccnductinc  toils  . 


Sine le- Sided  MMR 


Sing.e-sided  N'MR  systems  have  seen  developed  in  recent  years  to 
measure  moisture  Levels  in  structural  concrete  and  in  sucs-rface 
soils.  Like  Inside-Out  N’MR ,  discussed  aoeve,  tnese  systems  also 
utilize  two  of  t.he  components  discussed  in  Section  17.3.  First, 
an  electromagnet  wit.n  a  G-shaped  yoke  is  used  to  produce  tne  GC 
magnetic  field.  T.hen ,  a  spiral  rf  toil  is  used  for  ootn  trsiis- 
mitti.nc  t.ne  rf  oulse  and  receivi.ng  t.ne  GMR.  sianal.  A  s.<ett.n  of  a 
si.ngle-sided  system  of  this  type  is  shown  i.n  Figure  V- 1  2  ’  . 

The  Single-Sided  NMR  syste.m  for  moisture  detection  i.n  structural 
concrete  '3,4)  employs  an  iron-core  magnet  with  a  'J-snaped  yo)<e . 
The  interior  length  of  the  G-shaped  yoKe  is  11  inches,  tne  center 
to  ce.nter  distance  oetween  t.ne  magnet  pole  faces  is  accut  1: 
inc.nes  '38  cm),  and  the  syste.m  is  designed  to  operate  at  a  field 
strength  of  494  Gauss. 

The  useable  penetration  depth  of  this  syste.m  is  aoout  10  cm, 
although  the  magnet  produces  494  Gauss  fields  at  distances  of  u 
to  15  cm.  Presumably,  limitations  in  the  rf  components  preclud 
operation  at  penetration  depths  between  10  and  15  c.m .  .As  con¬ 
figured,  the  syste.m  operates  at  2.1  Megahertz  and  requires  2  00 
kilowatts  of  rf  power  to  produce  a  30  microsecond  90-degree 

p  U  1  S  0  • 


At  a  penetration  depth  of  2.5  cm,  the  sensitive  region  is 
approximately  7.5  cm  in  diameter  and  0.6  cm  thick.  As  expected, 
the  sensitive  volume  increases  with  increasing  distance  from  the 
instrument . 

Single-sided  syste.ms  designed  around  iron  magnets  wit.h  G-shaped 
yokes  are,  by  necessity,  quite  neavy .  The  e  lectrom.agnet  for  the 
system  described  above  weighs  approximately  600  pounds,  and  the 
associated  equipment  weighs  an  additional  450  pounds. 


C.  Advanced  Single-Sided  NMR  Syste.ms 

The  drawbacks  of  t.he  Inside-Out  MMR  geometry  for  single-siiec  GMR 
applications  are  oovious.  Furt.nermcre ,  it  would  be  difficult, 
per.naps  impossible,  to  engi.neer  a  s i.ng  1  e-s i ded  GMP.  system 
incorporating  a  G-shaped,  iron-core  m.agnet  cnac  wculc  exnioit  all 
tr. 0  d05irHci0  discussed  ir.  S^cticr.  A  cf  tr. is  ir. acesr. 
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Ma'cr  i.Tcrcvemer.Tis  to  present  tecnnclocy  are  required  in 
principal  areas  : 

1.  The  strenath  of  the  magnetic  field  in  tne  remote 
region  must  oe  increased. 

r.ne  penetration  depth,  or  distance  to  the  remote 
region  of  nomcgeneous  field,  must  be  increased. 

3.  Tne  coupling  of  the  rf  field  to  the  remote  region 
if  nomcgeneous  field  must  be  improved. 


To  address  tnese  issues,  three  potential  s i.ng  1  e- s ided  N.MR  con¬ 
figurations  were  considered.  The  first  utilizes  opposed  loops 
and  or  solenoids  to  produce  t.ne  re.mote  region  of  DC  field  and  a 
sem.itcroid  to  produce  t.he  rf  field.  T.his  is  termed  Opposed  Loop 
Se.mitoroid  .Re.mote  MMR . 

Tne  second  utilizes  a  se.mitoroid  to  produce  the  Ho  field  and 
oppcsed  loops  or  solenoids  to  produce  the  re.mote  rf  fielo.  This 
IS  termed  Se.mitoroid  Opposed  loop  Re.mote  MMR . 


The  tnird  utilizes  two  se.mitoroid 
tne  otner  for  the  rf  field.  This 
Re.mote  NMR .  These  t.hree  possible 
cussed  oelow. 


coils,  one  for  the  DC 
is  termed  Double  Se.mi 
syste.m  configurations 


oroid 
are  d i s - 


iDcosed  loop  Se.mitoroid  Sing  le-5ided  NMR  Svste.m 

.A.n  oppcsed  loop  se.mitoroid  si.ngle-sided  NMR  syste.m  is  shown  in 
Figure  7-4.  Steady  currents  flowing  t.hrough  the  oppcsed  loops 
produce  tne  re.mote  region  of  homogeneous  field.  Ho,  whicn  is 
directed  along  the  Z-axis.  The  rf  field.  Hi,  is  produced  by  a 
se.mitoroid  and  is  directed  perpe.ndicu  lar  to  the  Ho  field. 

An  advantage  of  t.nis  configuration  is  t.nat  tne  nomogenecus  region 
of  field  IS  large  and  well-defined.  However,  a  potential  draw- 
oacx  13  tnat  tne  rf  field  may  -.-ary  significantly  over  tne  -.’olume 
of  t.ne  Ho  field.  It  is  quite  .ixely,  t.nougn ,  tnat  field  snaping 
could  oe  ised  to  opt  i.mi  ze  coupling  oe  tween  t.ne  two  fields. 
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An  advantage  of  tnia  configuration  la  tnat  tne  rf  acil  *a  mcr- 
efficient  that  of  the  aemi toroid  opcosed  loop  ayste.r  auc’e. 
however,  tne  rf  coil  in  tnis  system  is  identical  tr  tnat  : :  tn 
opposed  loop  semitoroid  system,  and  tnat  acn f i lu r a t i an  naa  tne 
addit;.cnal  advantage  of  a  norocgeneous  DC  field. 


D.  Evaluation  of  Alternate  Remote  NMR  Syste.ms 

The  proposed  single-sided  NMR  systems  based  upon  opposed  ac,.ls  or 
solenoids  to  produce  the  DC  field  and  semitoroid  coils  to  pr 
tne  rf  field  possess  clear  qualititive  advantages  over  tne  o 
potential  configurations  discussed  in  the  previous  section.  In 
order  to  select  tne  most  attractive  of  the  re.maining  systems,  tne 
f  igure-of -merit ,  introduced  in  Cnapter  17,  was  e.mplcyed. 

The  f igure-of -merit ,  ?s ,  was  calculated  for  eacn  of  tne  ro^r 
possiole  systems.  Since  se.mitoroids  are  itiliced  as  rf  coils, 
t.ne  systems  differ  i.n  t.he  met.hods  used  to  produce  tne  EC  fields. 
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calculated  value  of  Fs  is  58.9  x  '£-3'.  As  the  sffsst  is 
increased  to  S  =  +2.5,  Fs  falls  to  6.56  x  (E-3),  almost  a  factor 
of  ten. 

However,  whe.n  the  solenoids  are  offset  in  t.he  opposite  direction 
to  S  =  -2.5,  Fs  increases  to  7.7  x  (E-5).  This  represents  an 
i.ncrease  of  over  a  factor  of  100.  T.he  latter  value  of  Fs 
represents  a  significant  improve.ment  over  any  of  the  previous 
configurations . 


Inside-Out  He  Imho 1 tz/Semi toroid  System 


Calculations  of  Fs  and  ot.her  relevant  parameters  for  Inside-Out 
Hel.mholtz/ Se.mitoroid  systems  are  presented  in  Table  V-4.  Al- 
t.nough  the  maximum  attainable  field  strengths  for  these  syste.ms 
are  comparable  to  those  of  the  previous  configurations,  tne 
re.mote  regions  of  homogeneous  field  are  much  larger.  This  leads 
to  a  direct  increase  in  the  f igure-of-meri t . 

For  example,  with  an  inner  coil  diameter  of  five,  a  set  of  simple 
opposed  loops  gives  Fs  =  53.3  x  (E-8).  The  comparable  value  for 
an  Inside-Out  Helmholtz  configuration  is  Fs  =  22.5  x  {E-5),  whicn 
IS  over  a  factor  of  400  larger.  In  addition,  this  Fs  is  al.most  a 
factor  of  three  larger  than  the  most  favorable  configuration  of 
offset  opposed  solenoids. 

The  f igure-of-merit  is  a  strong  function  of  the  diameter  of  the 
inner  coil,  and  as  this  diameter  increases,  Fs  decreases.  How¬ 
ever,  both  the  overall  size  of  the  coil  system  and  the 
penetration  depth  increase  with  increasing  inner  coil  diameter, 
so  practical  considerations  may  lead  to  compromises  in  the  choice 
of  syste.m  configuration. 


Taole  7-4 

f:gl’re-of-merit  and  characteristics 
OF  inside-out  hel:^.holtz  opposed  coils 
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The  point  of  maxi.mum  field  on  the  Z  axis  is  Zm .  The  field  strengtn 
at  Z.m  13  3m.  The  diameter  of  the  homogeneous  region,  defined  by  a 
100  ppm  cna.nge  in  3,  is  dZ  and  the  volume  of  this  region  is  Vm. .  Mo 
IS  tne  nuclear  magnetization,  and  Fs  is  the  f igure-of -men t . 
Notations  of  t.ne  form  x'E  +  Y)  mean  "ti.mes  ten  raised  to  the  Y-th 
power  "  . 
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A.  Introducrion 

The  f igure-of-merit  evaluations  of  tne  previous  chapter  lead 
directly  to  a  preferred  configuration  for  an  advanced  single¬ 
sided  MMR  system  for  remote  detection  appl ications .  Specific 
features  of  the  preferred  configuration,  estimated  signal -to- 
noise  ratios,  and  practical  considerations  associated  v/ith 
implementation  are  discussed  in  this  chapter. 


3.  Syste.m  Configuration 

Based  upon  the  f igure-of -merit  calculations,  it  appears  that  a 
single-sided  remote  NMR  system  t.hat  utilizes  Inside-Out  Hel.mholtz 
coils  to  produce  the  steady  Ho  field  possesses  clear  advantages 
over  competing  configurations.  With  proper  design,  ‘suc.n  a  system 
snould  exhibit  all  or  most  of  the  desiraole  attributes  discussed 
in  the  previous  chapter,  i.e.,  a  reasonable  penetration  depth, 
portaoility  or  transportability,  and  a  well-defined  remote  region 
of  inspection. 

A  sketch  of  the  proposed  system  is  shown  in  Figure  VI-1.  In  tnis 
sketch,  the  semi toroid  is  shown  with  a  major  diameter  approxi¬ 
mately  equal  to  the  diameter  of  the  inner  coil  of  the  Inside-Cut 
Helmholtz  pair.  The  optimum  configuration  would,  of  course, 
depend  upon  the  particular  choice  of  Inside-Out  Hel.mholtz 
geometry.  However,  assuming  Ra  =  2 ( Rb ) ,  the  configuration 
illustrated  would  provide  an  rf  field  strengt.h  at  the  ce.nter  of 
the  re.mote  region  approximately  20  -  25%  of  the  field  strengt.h 
available  at  t.he  semitoroid  face.  In  addition,  the  inhomogenei ty 
of  the  rf  field  over  the  remote  volume  is  estimated  to  be  less 
than  10%,  which  is  very  attractive. 

It  may  be  desirable  to  design  the  system  with  two  separate  rf 
coils,  one  for  transmitting  and  the  other  for  receiving.  Such  a 
design  may  lead  to  a  more  efficient  use  of  transmitter  power  or 
to  an  i.mprovement  in  the  signal-to-noise  ratio  t.hrough  reduced 
system  recovery  time,  improved  detector  coil  Q,  or  the  like. 

If  separate  receiver  and  transmitter  coils  are  required,  a  si.mple 
solution  would  be  to  incorporate  a  solenoid  detector  coil  as 
shown  i.n  Figure  71-2.  The  detector  coil  would  t.hen  be  oriented 
to  minimize  coupling  witn  t.he  semi  toroid.  An  attractive  alter¬ 
native  may  be  to  e.mploy  crossed  semitoroids  si.milar  to  those 
shown  in  Figure  V-6  (24). 
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C.  £  igna  I -Tc-N'oise  Ratio 

Equation  II-"  was  used  to  calculate  tne  s  i  ma  x  - 1  c-r.c  i  se  ratio 
expected  to  oe  observed  from  a  sincle-sided  N'>!R  system  prcducing 
a  1,300  Gauss  field  in  a  remote  region  10  centimeters  away. 

Under  the  following  conditions. 


Proton  Density  =  6.6  x  10^^; 
Sample  Temperature  =  300  Kelvin; 
System  Bandwidth,  B  =  42.57  KHz; 
Ra  =  2 ( Rb )  =  40  centimeters 
d  =  Zm ;  Z  =  3Zm/2;  ,^/d  =  0.5; 


and  with  the  factor  F  set  equal  to  unity,  the  SNR  was  calculated 
to  be 


SNR  =  1 1  .-4  . 


Although  there  are  a  number  of  effects  that  will  act  to  reduce 
the  expected  SNR  from  the  calculated  value,  this  is  an  encourag¬ 
ing  result.  In  addition,  it  is  likely  that  the  above  conditions 
could  be  optimized  to  produce  an  even  greater  estimated  SNR. 

An  interesting  result,  which  arises  from  tne  geometry  of  the 
Inside-Out  Helmholtz  system,  is  that  the  estimated  SNR  for  a 
particular  configuration  is  independent  of  the  penetration  depth. 
That  is,  if  the  size  of  the  above  system  is  doubled,  and  the  Ho 
field  in  the  remote  region  remains  the  same,  then  the  predicted 
SNR  is  again  11.4.  The  reason  is  t.nat  the  volume  of  the  remote 
region  varies  as  tne  cube  of  the  penetration  depth.  T.herefore, 
at  constant  Ho  field,  the  signal -to-noise  ratio  is  independent  of 
penetration  depth.  Of  course,  practical  constraints  on  coil 
size,  rf  field  strength  and  the  like,  make  it  very  unlikely  that 
the  SNR  would  be  penetration  depth  independent,  and  these  ef feces 
are  addressed  below. 

It  snould  also  be  noted  that  ehere  are  effects  whicn  will  act  to 
reduce  the  signal -to-noise  ratio  calculated  above.  First,  the 
SNR  depends  directly  upon  the  proton  density,  and  a  .material  witn 
a  high  proton  density  was  assumed.  Hov/ever,  mar.v  e.xplosive 
compounds  also  have  high  proton  densities,  so  this  may  not  pose  a 
serious  oroble.m. 

Second,  tne  SNR  calculated  above  assumes  that  F  =  1.3.  The 
quantity  F  is  a  factor  that  i.ncludes  the  effects  of  finite  rf 
pulse  length,  rf  field  inhomogeneity,  syste.m  recovery  ti.me  and 
signal  loss  due  to  short  T2  times.  Small  changes  i.n  rf  pulse 


ler.ctr.  or  sysrem  recovery  tine  ter.  leaa  to  t or  r. - 
i.t  tne  ODserved  SNR.  Si.nce  nost  e.xpicsive  cortpcunds  ex.ticit 
snort  T2  ti.nes,  tnese  effects  are  likely  to  oe  siynifioant. 

r 


T  a(  Z  ,  ror  Z  ■  <  Zo; 

F  tf<  Z  for  Z  > ->  Zc; 


It  would  be  difficult  to  deterni.ne  the  factor  ?  quantitati'/e 
out  qualitative  esti.tates  are  fairly  straightforward.  F  is 
lixely  to  oe  t  ■..-ax  1  y- varyi.nc  fu.ncticn  of  pe.netrat ion  dept;'. 
s~all  Z  and  a  s tronc  1  y-vary mg  function  of  penetration  dept.i 
large  Z.  A  reasonable  as  su.t.ption ,  based  upon  e.xperience  wit 
I.nside-Out  N'MR  syste.n.s  iTjZS,!,  is  t.nat 


wnere  Zo  is  a  characteristic  distance  dencti.ng  t.ne  cna.nce  i.n 
regions . 


An  ana...ysis  or  tne  electrical  and  pnysical  characteristics  of  the 


quire.nents  can  be  .net  for  coil  diameters  between  about  1  and  2  7 
centi.meters .  While  the  lower  li.mit  does  not  present  a  serious 
constraint,  the  upper  limit  suggests  that  the  transition  region 
occurs  near  Zo  =  27  cm. 


The  syste.m  recovery  time  depe.nds,  to  some  extent,  upon  the  Larmor 
frequency,  or  static  field,  chosen  for  observation.  The  higher 
t.ne  frequency  or  field,  in  general,  the  shorter  the  recovery 
ti.T.e .  In  addition,  the  shorter  the  T2  of  the  material  under 
study,  t.ne  greater  the  loss  of  signal  during  the  recovery  time. 

These  considerations  were  collected  into  the  following  ad  hoc  e.x- 
pression  for  the  factor  F, 


=  D 


exp{-l/Ho)/( iZ/Zo) 


2/Zo)  ) , 


'vnere  D  and  Zo  are  constants  and  Ho  is  the  applied  field  in 
xilogauss.  Since  th  T2  times  of  most  explosive  compounds  are 
si.milar,  and  short,  the  dependence  of  F  on  T2  was  not  explicitly 
included.  I.mplicit  in  this  expression  is  the  assumpticn  that  tne 
effect  of  t.ne  recovery  time  will  reduce  tne  available  5XR  oy  a 
factor  of  'l/e'i  at  a  field  of  1,000  Gauss,  and  more  at  lower 
fields  . 


(t  (t 


The  effects  ef  rf  field  inncmcceneity  and  ccapling  inefficiencies 
were  tnen  assumed  to  further  reduce  the  SNR  by  a  factor  of  two  at 
a  field  of  1,1)00  Gauss  and  a  penetration  deptn  of  10  cm.  The 
factor  ?  was  tnen  evaluated  under  the  conservative  assumption,  So 

—  '0  ""Ti 

Wit.n  tne  inclusion  of  these  numerical  val..es,  £quat_or.  '/I-l 
oecomes : 


R  =  =  S  1.0 
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wnere  agai.n  Ho  is  the  applied  field  in  ;<ilogauss. 

Equations  II-7  and  "71-2  were  used  to  estimate  s igna  1  - to-noise 
ratios  for  a  variety  of  penetration  depths  and  applied  fields. 

The  results  are  presented  in  Table  VI-1,  and  the  strong  depen¬ 
dence  of  SNR  on  applied  field  is  immediately  obvious.  At  0.5 
hilogauss,  a  penetration  depth  of  2  cm  provides  a  barely  adequate 
SNR,  while  at  1.5  kilogauss,  penetration  depths  of  50  cm  appear 
feasible . 

The  SNR  also  depends  strongly  upon  the  configuration  of  the 
I.nside-Out  Hel.mholtz  coils.  Predicted  signal -to-noise  ratios  for 
a  variety  of  configurations  are  presented  in  Table  VI-2.  Two 
cases  are  covered  in  this  table.  In  the  first,  the  Ho  fields  in 
the  re.mote  regions  are  assumed  to  vary  as  indicated  in  Table  V-4, 
with  a  1,000  Gauss  field  for  the  Ra  =  40  cm  system.  In  the 
second,  a  1,000  Gauss  field  is  assumed  for  all  configurations. 

It  IS  obvious  from  Table  VI-2  that  the  SNR  is  strongly  dependent 
upon  t.ne  syste.m  configuration.  Unfortunately,  t.he  most  attrac- 
ive  configurations,  from  a  practical  standpoint,  also  exhibit 
he  lowest  signal-to-noise  ratios. 
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SINGLE-SIDED  NMR  SIGNA: 
VERSUS  PENETRATION  DEPTH 


,-TO->iCISE  R.' 
AND  APPLIED 


Signal -to-noise  rarios  for  aincle-sided  remote  detection 
as  a  function  of  the  penetration  depth  in  centimeters  and 
field  in  kilogauss.  An  Inside-Out  Helmholtz  geometry  wit 
IS  assumed. 


SIGNAL-TO-NCISE  RATIO 


m  ( cm ) 

0  .  5  KG 

1.0  KG 

1.5  KG 

2 

0.497 

5.40 

17.0 

5 

0  .  304 

3 .30 

10.4 

10 

0.  193 

2.09 

6.58 

20 

0 .096 

1.05 

3.29 

30 

0.053 

0.575 

1.81 

40 

0.031 

0.341 

1  .07 

50 

0 . 020 

0.217 

0 . 682 

60 

0.014 

0.  147 

0 . 461 

NMR  sy 
tne  s 
h  Ra  = 


^  tl  ( 1 
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Signal -to-ncise  ratios  in  sincle-sided  systems  emtloyi.ty  Ir.site- 

Cat  Helmncitt  coiis.  As  indicated  below,  one  penetration  aeptr.  is 
centimeters,  and  tne  diameter  of  the  inner  coil  ^  Rit  '  is  20  centim.e- 
ters.  The  diameter  of  the  outer  coil  varies  depending  upon  tne 
particular  Inside-Out  Helmholtz  configuration  chosen. 

3i gna 1 - tc-noise  ratios  are  presented  for  two  separate  rases,  A  ana 
3;.  The  3 igna  1  - to-noise  ratios  presented  in  colum.n  five.  Case  A, 
assume  t.nat  the  a.mpere- turns  in  tne  large  diameter  outer  coil  is  nelu 
constant,  i.ndependent  of  configuration,  ana  t.nat  the  field  stre.nct:!  in 
the  re.mcte  region  varies  as  indicated  in  Taoie  V-4.  Furtnermore,  for 
Ra  =  40  centi.meters  ,  the  static  field  is  assumed  to  oe  1,000  Gauss. 

The  3  igna  1  -  to-noise  ratios  presented  i.n  column  six.  Case  3,  assu.me 
t.ne  field  stre.ngth  i.n  the  remote  region  is  i.ndepe.nde.nt  of 
iguration.  A  field  of  1,000  Gauss  is  assu.med  for  eacn  rase. 


Ra 

Rb 

Zm 

Ho 

SNR  (A) 

SNR  ■  3 

66.' 

20  .  0 

10.0 

1 .40 

10.5 

4.02 

50  .  0 

20 , 0 

10.0 

1.30 

4.74 

2  .  ~9 

40.0 

20  .  0 

10.0 

1  .00 

2.09 

2.09 

3  3  * 

20 . 0 

10.0 

0.799 

0.309 

1.63 

23.6 

20 . 0 

10.0 

0.601 

0.242 

i.  .  J  1 

2  5.1 

20.0 

10.0 

0.399 

0.038 

1.08 

. r.  zc  irr'.r: . -irr.er. ir.e  r;rccos6G  s  ir. -  .  s i  ^ifrc  ...■': 

r.jrr.cer  c:  ^rac*ical  consi  ieratior.s  t.zsz  ce  adcrass' 
particular,  tne  currer.ts  required  to  produce  one  c 
tne  pr.ysical  ad'/antaces  and  1  i.Titations  or  specif  i 
uecrr(etries ,  and  any  electrical  electronic  li.T.itati, 
t a /< e n  into  account.  T .n e s e  topics  are  a i s c u s 5 e c  in 
sections  union  follow  oelcw. 


:urr' 


.Recuire.tients 


T.te  requirea  numcer  or  ampere-turns  to  proauce  a  gz'/er.  :ie-a 
stre.ngt.n  i.n  t.ae  re.mote  regio.n  was  calculated  for  several  Inside- 
Out  -He  l.m.no  1 1 2  coil  geometries.  The  results  for  one  particular 
geometry,  i.e.  ,  the  .Ra  =  2  { Rb )  geometry,  are  shown  in  Figure 
VI  -3.  This  13  a  plot  of  t.he  ampere-turns  required  in  the  large 
diameter  coil  as  a  function  of  t.he  penetration  depth,  2m,  and  tne 
field  stre.ngt.n  in  the  remote  region.  Ho.  It  can  oe  seen  t.nat 
well  over  100,000  a.mpere- turns  are  required  to  produce  a  re.mote 
field  of  1,000  Gauss  at  a  distance  of  10  cm. 


The  required  number  of  ampere-turns  depends  upon  the  specific 
geo.metry  of  the  Inside-Out  Helmholtz  pair.  For  constant 
penetration  depth  and  remote  region  field  strength,  t.he  requir 
number  of  ampere-turns  in  the  outer  coil  is  a  minimum  for  the 
=  geometry.  This  is  s.hcwn  in  Figure  Vl-<!,  unhich  assumes 

outer  coil  diameter  Ra  =  10  units. 


The  required  number  of  ampere-turns  for  the  inner  coil  of  tne 
Inside-Cut  Helmholtz  pair  is  shown  in  Figure  VI-5.  The  calcu¬ 
lations  for  this  plot  also  assu.me  an  outer  coil  diameter  of  10. 

From  tne  results  shown  in  Figures  VT-3,  VI-4  and  VI-5,  it  is 
obvious  t.nat  it  would  be  very  difficult,  and  likely  i.m.possio  le , 
to  produce  substantial  fields  at  significant  pe.netration  depths 
using  conventional  resistive  coil  technology.  Pulsed  fields 
could  oe  used  in  certain  applications,  but  t.he  resulting  syste.m 
would  li.kely  be  quite  complex,  especially  since  explosives  tend 
to  exhibit  long  T1  times. 

E.mployi.ng  super conducti.ng  coils  .Tiay  oe  an  attractive  sc.uticn  nc 
this  prcole.m  since  it  is  li.iely  t.nat  super  conduct  i..nc  rcils  iculc 
produce  iseaole  re.mote  f.eld  strengths  and  penetration  depths. 


course,  t.nere  are  drawcac.-cs  to  t.ne  use  of  sucerconcuctinc 
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Cell  Geometries 


The  Inside-Out  Helmholtz  configuration  offers  a  ide  range  of 
geometries  suitable  for  single-sided  NMR.  This  range  is 
illustrated  in  the  sketches  of  Figures  VI-6  through  VI-11. 

These  sketches,  which  are  drawn  to  scale,  show  six  separate 
Inside-Out  Helmholtz  geometries.  The  diameter  of  the  outer  coil 
IS  held  constant  in  each  of  the  sketches.  The  spheres  on  the 
Z-axis  represent  the  relative  sizes  of  the  remote  regions  of 
field,  with  a  1%  change  in  field  representing  the  boundary. 

It  IS  obvious  from  this  series  of  sketches  that  some  coil 
geometries  are  more  attractive  than  others  for  single-sided  MMR. 
As  the  diameter  of  the  inner  coil  increases,  the  distance  to  tne 
re.mote  region  (i.e.,  the  penetration  depth)  also  increases.  In 
addition,  the  size  of  the  remote  region  increases.  Therefore, 
those  geometries  with  relatively  large  inner  coil  diameters  would 
appear  to  be  very  attractive  for  single-sided  MMR  applications. 

However,  as  the  inner  coil  diameter  increases,  the  figure-of- 
merit  for  the  system  decreases.  This  is  a  real  effect  and  means 
t.nat,  all  things  being  equal,  a  system  with  a  smaller  inner  coil 
diameter  will  outperform  one  with  a  larger  inner  coil  diameter. 

There  are  two  reasons  for  this  result.  First,  the  coupling 
between  the  remote  region  and  the  detector  coil  (and,  hence,  the 
f igure-of-merit )  varies  as  the  inverse  cube  of  the  penetration 
depth.  Second,  the  strength  of  the  Ho  field  in  the  remote  region 
decreases  with  increasing  penetration  depth. 

The  first  effect  is  common  to  all  re.mote  MMR  configurations  and 
little  can  be  done.  Since  there  may  be  little  advantage  in 
operating  the  system  at  fields  above  about  1,000  Gauss,  because 
of  level  crossing  effects,  the  second  effect  may  be  partially 
compensated  by  choosing  an  attractive  geometery  and  then 
increasing  the  current  in  the  coils  to  produce  the  required 
external  field.  However,  this  strategy  cannot  be  pursued 
indefinitely  since  the  coils  have  current  li.mits.  These  limits 
are  discussed  in  the  following  section. 


Coil  Li.mitations 


There  are  two  potentially  serious  drawbacks  to  the  use  of  super¬ 
conducting  coils  in  t.he  proposed  system.  First,  superconducting 
coils  .must  be  immersed  in  liquid  helium.  This  problem  would  be 
very  difficult  to  overcome  were  it  not  for  a  recent  advance  in 
cryogenic  tec.hnology:  t.he  development  of  insulated  containers 
that  require  refilling  only  three  or  four  times  a  year  and  liccle 
.mainte.nance .  For  a  stationary  field  instrument,  this  should  be 
adequate.  For  a  portable  system,  thougn ,  the  requirement  for 
liquid  helium  may  still  pose  serious  problems. 
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Second,  the  maximum  remote  field  attainaoie  from  a  superccna^c- 
ting  Inside-Out  Helmnoltz  pair  may  not  oe  adequate.  The  maxi.mum 
field  depends  upon  the  maximum  current  chat  can  flow  througn  t.ce 
coils,  and  the  maximum  current  is  determined  oy  tne  magnitudes  of 
the  fields  produced  at  the  sites  of  the  coil  windings. 

Superconduct i.ng  materials  can  only  tolerate  magnetic  fields  up  cc 
a  certain  .maximum  before  reverting  to  their  normal  state.  The 
fields  at  the  sites  of  the  coils  of  an  Inside-Out  .Hel.mholtz  pair 
are  determined  by  the  configuration  chosen  and  by  the  curre.nts  in 
the  coils.  The  maximum  remote  field  is  reached  wnen  the  li.miting 
current  is  reached  in  one  or  both  of  the  coils. 


To  address  this  problem,  we  contacted  Dr.  James  Carolan,  Pres¬ 
ident  cf  the  Nalorac  Cryogenics  Corporation,  a  manufacturer  of 
superconducting  coils  and  associated  components.  Dr.  Carolan 
analyzed  this  proble.m  using  a  computer  program  developed  for  this 
purpose  (26).  His  conclusions,  for  the  Sa  =  2 ( Rb )  geometry,  were 
that  a  1,00G  Gauss  remote  field  with  a  penetration  depth  of  30  cm 
IS  certai.nly  possible  and  that  remote  fields  up  to  2,000  Gauss 
m.av  oe  achievaole. 
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VII.  CONCLUSIONS  AND  RECOMMENDATIONS 


A.  Conclusions 

The  results  of  this  study  lead  to  the  following  principal 
conclusions  regarding  re.T’.ote  detection  of  explosives  using 
single-sided  NMR  technology. 


1.  Re.mote  detection  by  single-sided  NMR  appears 
feasible.  The  present  analysis  indicates  tnat 
useable  penetration  depths  and  adequate  signal- 
to-noise  ratios  are  attainable. 

2.  NMR  remote  detection  of  certain  explosives  is 
possible,  but  detailed  data  on  the  NMR  proper¬ 
ties  of  explosives  are  sparse.  It  may  also  be 
possible  to  define  a  unique  signature,  using 
hydrogen/nitrogen  level  crossing,  for  each 
explosive  compound.  Additional  data  and 
analyses  are  required. 

3.  An  advanced  NMR  remote  detection  system  will 
require  superconducting  coils  to  produce  the 
.magnetic  field. 

4.  The  preferred  configuration  will  employ  an 
Inside-Out  Helmholtz  pair  to  produce  the 
magnetic  field,  and  a  semitoroid  coil  to 
produce  the  rf  field. 


B.  Recommendations 

T.hese  conclusions  lead  directly  to  the  two  recommendations  pre¬ 
sented  below. 


Recommendation  1 .  Acquire  NMR  data  on  the  compounds 
of  direct  i.nterest.  Proton  relaxation  times  ana 
the  effects  of  hydrogen/nitrogen  interactions  are 
of  particular  importance. 


Recommendation  2 .  Perform  a  feasibi 
f leld-useable  system.  Use  t.he 
ation  1,  above,  to  address  the 
discrimination  between  various 
i.nspection  times  u.nder  various 


lity  study  for  a 
data  from  Reconmend- 
possibility  of 
explosives.  Esti.mate 
sets  of  conditions. 
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VIII.  NOMENCLATURE 


Correlation  Ti.nie.  For  random  processes  such  as  molecular 

m.cticn,  a  measure  of  the  time  dependence  of  the  fluctuations. 

In  NMR,  correlation  times  near  the  inverse  of  the  Larmor 
frequency  produce  short  T1  times. 

Dead  Ti.me.  The  ti.me  for  the  NMR  system  to  recover  followi.nc 
an  rf  pulse,  i.e.,  the  time  it  takes  for  an  FID  to  be  visiole 
after  a  pulse.  Also  called  the  recovery  time. 

FID .  The  NMR  free  induction  decay  signal  following  t.he 

application  of  an  rf  pulse,  or  pulses,  at  the  Larmor  frequency. 

Free  Induction  Decay.  See  FID  above. 

Ho  Field.  The  externally-applied  magnetic  field  responsible 
for  nuclear  precession  at  the  Larmor  frequency. 

HI  Field.  The  rf  field  applied  at  the  Larmor  frequency  to  tip 
the  nuclear  magnetization.  The  HI  field  is  applied  perpendi¬ 
cular  to  t.he  Ho  field. 

Inside-Out  Helmholtz.  A  set  of  coils  with  spacing,  radii  and 
current  flows  arranged  to  produce  a  remote  region  of  homoge¬ 
neous  field.  At  the  remote  field  maximum,  both  the  first  and 
second  derivatives  of  field  strength  with  respect  to  distance 
are  zero. 

Inside-Out  NMR.  A  particular  type  of  remote  NMR  configuration 
in  which  opposed  magnets  are  used  to  produce  a  toroidal  remote 
region  of  homogeneous  field.  Originally  developed  for  well 
logging  applications. 

Larmor  frequency.  The  frequency  at  which  the  NMR  experi.ment  is 
performed;  the  precession  frequency  of  the  nuclei  in  the  static 
magnetic  field.  The  Larmor  frequency  is  proportional  to  the 
strengt-h  of  the  static  magnetic  field  and,  for  protons,  is  4.26 
.MHz  at  1,000  Gauss. 

Polarization  Time.  The  duration  of  time  that  the  item  to  be 
inspected  is  immersed  i.n  t.he  static  magnetic  field  before 
interrogation  by  the  NMR  pulse  sequence. 

Recovery  Time.  See  dead  time  above. 

RF .  Radio- : requency ;  in  NMR,  frequencies  in  the  range  of 
about  0.3  -  3*00  .MHz. 

RF  Pulse.  A  short,  inte.nse  burst  of  rf  energy  applied  tc  the 
nuclei  by  the  transmitter. 
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RF  Pulse  Sequence.  A  sequence  of  rf  pulses  designed  to  elicit 
a  specific  response  from  the  nuclear  spin  system. 

Solid  Echo.  A  response  from  stongly  coupled  nuclei  in  solids 
produced  oy  a  sequence  of  two  closely  spaced  rf  pulses.  The 
solid  echo  is  often  used  to  improve  the  SNR  in  syste.ms  while 
observing  signal',  from  short  T2  materials. 

SNR.  T.he  signal-to-noise  ratio  of  the  NMR  system;  usually 
refers  to  the  voltage  SNR. 

Soin-Lattice  Relaxation  Time.  See  Tl. 


Soin-Soin  Relaxation  Time.  See  T2 . 

Static  Field.  See  Ho  field  above. 

Tl .  The  spin-lattice  relaxation  time  of  the  nuclear  spin 

syste.m.  The  characteristic  time  for  the  nuclear  spin  system 
to  return  to  thermal  equilibrium  with  its  surroundings 
following  a  disturbance. 

T2 .  T.he  spin-spin  relaxation  ti.me  of  the  nuclear  spin  syste.m. 
The  characteristic  ti.me  for  the  nuclear  spin  syste.m  to  rerurn 
to  internal  equilibrium  following  a  disturbance. 
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